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Solar  radiation  is  a  primary  driver  for  many  physical,  chemical,  and  biological  processes  on  the  earth’s 
surface.  Solar  energy  engineers,  architects,  agriculturists,  hydrologists,  etc.  often  require  a  reasonably 
accurate  knowledge  of  the  availability  of  the  solar  resource  for  their  relevant  applications  at  their  local. 
In  solar  applications,  one  of  the  most  important  parameters  needed  is  the  long-term  average  daily 
global  irradiation.  For  regions  where  no  actual  measured  values  are  available,  a  common  practice  is  to 
estimate  average  daily  global  solar  radiation  using  appropriate  empirical  correlations  based  on  the 
measured  relevant  data  at  those  locations.  These  correlations  estimate  the  values  of  global  solar 
radiation  for  a  region  of  interest  from  more  readily  available  meteorological,  climatological,  and 
geographical  parameters.  The  main  objective  of  this  study  is  to  chronologically  collect  and  review  the 
extensive  global  solar  radiation  models  available  in  the  literature  and  to  classify  them  into  four 
categories,  i.e„  sunshine-based,  cloud-based,  temperature-based,  and  other  meteorological  parameter- 
based  models,  based  on  the  employed  meteorological  parameters  as  model  input. 

Furthermore,  in  order  to  evaluate  the  accuracy  and  applicability  of  the  models  reported  in  this  paper 
for  computing  the  monthly  average  daily  global  solar  radiation  on  a  horizontal  surface,  the  geographical 
and  meteorological  data  of  Yazd  city,  Iran  was  used.  The  developed  models  were  then  evaluated  and 
compared  on  the  basis  of  statistical  error  indices  and  the  most  accurate  model  was  chosen  in  each 
category.  Results  revealed  that  all  the  proposed  correlations  have  a  good  estimation  of  the  monthly 
average  daily  global  solar  radiation  on  a  horizontal  surface  in  Yazd  city,  however,  the  El-Metwally 
sunshine-based  model  predicts  the  monthly  averaged  global  solar  radiation  with  a  higher  accuracy. 
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1.  Introduction 

Solar  radiation  arriving  on  earth  is  the  most  fundamental 
renewable  energy  source  in  nature.  Solar  energy,  radiant  light 
and  heat  from  the  sun,  has  been  harnessed  by  humans  since 
ancient  times  using  a  range  of  ever  evolving  technologies.  The 
energy  from  the  sun  could  play  a  key  role  in  de-carbonizing  the 
global  economy  alongside  improvements  in  energy  efficiency  and 
imposing  costs  on  greenhouse  gas  emitters.  In  the  studies  of  solar 
energy,  data  on  solar  radiation  and  its  components  at  a  given 
location  are  a  fundamental  input  for  solar  energy  applications 
such  as  photovoltaic,  solar-thermal  systems,  solar  furnaces  and 
passive  solar  design.  The  data  should  be  reliable  and  readily 
available  for  design,  optimization  and  performance  evaluation  of 
solar  systems  for  any  particular  location.  The  best  way  to 
determine  the  amount  of  global  solar  radiation  at  any  site  is  to 
install  measuring  instruments  such  as  pyranometer  and  pyrheli- 
ometer  at  that  particular  place  and  to  monitor  and  record  its 
day-to-day  recording,  which  is  really  a  very  tedious  and  costly 
exercise  [1],  In  spite  of  the  importance  of  solar  radiation  mea¬ 
surements,  in  many  developing  countries  this  information  is  not 
readily  available  because  of  not  being  able  to  afford  the  measur¬ 
ing  equipments  and  techniques  involved.  Therefore,  a  number  of 
correlations  and  methods  have  been  developed  to  estimate  daily 
or  monthly  global  solar  radiation  on  the  basis  of  the  more  readily 
available  meteorological  data  at  a  majority  of  weather  stations. 
Empirical  models  which  have  been  used  to  calculate  solar  radia¬ 
tion  are  usually  based  on  the  following  factors  [2]: 


(1)  Astronomical  factors  (solar  constant,  earth-sun  distance,  solar 
declination  and  hour  angle). 

(2)  Geographical  factors  (latitude,  longitude  and  elevation  of  the 
site). 

(3)  Geometrical  factors  (azimuth  angle  of  the  surface,  tilt  angle  of 
the  surface,  sun  elevation  angle,  sun  azimuth  angle). 

(4)  Physical  factors  (scattering  of  air  molecules,  water  vapor  content, 
scattering  of  dust  and  other  atmospheric  constituents  such  as  02, 
N2,  C02,  O,  etc.). 

(5)  Meteorological  factors  (extraterrestrial  solar  radiation,  sun¬ 
shine  duration,  temperature,  precipitation,  relative  humidity, 
effects  of  cloudiness,  soil  temperature,  evaporation,  reflection 
of  the  environs,  etc.). 

The  number  of  correlation  that  have  been  published  and  tested 
to  estimate  the  monthly  average  daily  global  solar  radiation  is 
relatively  high,  which  makes  it  difficult  to  choose  the  most 
appropriate  method  for  a  particular  purpose  and  site.  Selecting 
an  appropriate  methods  from  various  existing  models  is  based  on 
their  data  requirements  (the  selected  methods  utilize  only  daily 
variables  normally  available  at  a  majority  of  weather  stations) 
and  the  model  accuracy. 

Estimations  of  the  monthly  mean  daily  global  solar  radiation 
for  a  large  number  of  locations  are  presented  in  various  research 
works.  The  main  objective  of  this  study  is  to  comprehensively 
collect  and  present  the  global  solar  radiation  models  available  in 
the  literature,  including  the  study  carried  out  on  the  estimation  of 
the  monthly  average  daily  global  solar  radiation  on  horizontal 
surfaces  and  to  classify  them  into  four  categories  based  on  the 


employed  meteorological  parameters  as  model  inputs.  It  is 
noteworthy  that  all  the  proposed  models  contain  empirical 
constants  which  depend  on  the  season  and  the  geographical 
location  of  a  particular  place.  The  empirical  coefficients  for  some 
correlations  available  in  the  literature  are  presented.  The  models 
collected  and  reviewed  in  this  research  are  chronologically  pre¬ 
sented  and  therefore  are  useful  for  selecting  the  appropriate 
model  to  estimate  global  solar  radiation  for  a  particular  place  of 
interest. 

In  addition,  in  order  to  evaluate  the  performance  of  the  models 
in  each  of  the  four  categories  described  in  Section  2,  some 
selected  models  from  each  category  are  used  to  estimate  the 
monthly  average  daily  global  solar  radiation  by  using,  astronom¬ 
ical,  geographical,  geometrical  and  meteorological  data  of  Yazd 
city,  Iran.  The  solar  radiation  values  produced  by  the  selected 
models  are  then  compared  with  each  other  as  well  as  with  the 
available  experimental  data.  Among  the  selected  models  in  each 
category,  the  model  which  produces  the  most  accurate  values  of 
global  solar  radiation  is  selected  and  used  for  predicting  the 
monthly  average  daily  global  solar  radiation  in  Yazd. 

2.  Modeling  of  global  solar  radiation 

Solar  researchers  have  developed  many  empirical  correlations 
which  determine  the  relation  between  solar  radiation  and  various 
meteorological  parameters.  As  the  availability  of  meteorological 
parameters,  which  are  used  as  the  input  of  radiation  models  is  the 
most  important  key  to  choose  the  proper  radiation  models  at  any 
location,  empirical  models  can  be  mainly  classified  into  four  following 
categories  based  on  the  employed  meteorological  parameters: 


(1)  Sunshine-based  models. 

(2)  Cloud-based  models. 

(3)  Temperature-based  models. 

(4)  Other  meteorological  parameter-based  models. 

Among  all  such  meteorological  parameters,  bright  sunshine 
hours,  cloud  cover  and  temperature  are  the  most  widely  and 
commonly  used  ones  to  predict  global  solar  radiation  and  its 
components  at  any  location  of  interest. 

2.1.  Sunshine-based  models 

The  most  commonly  used  parameter  for  estimating  global 
solar  radiation  is  sunshine  duration.  Sunshine  duration  can  be 
easily  and  reliably  measured,  and  data  are  widely  available  at  the 
weather  stations.  Most  of  the  models  for  estimating  solar  radia¬ 
tion  that  appear  in  the  literature  only  use  sunshine  ratio  (S/S0)  for 
prediction  of  monthly  average  daily  global  radiation.  In  the 
following  section,  correlations  which  use  only  the  sunshine  ratio 
as  the  key  input  parameter  are  presented  and  classified  based  on 
their  developing  year. 

2.1.1.  Model  1 :  Angstrom-Prescott  model 

The  first  and  the  most  widely  used  correlation  for  estimating 
monthly  average  daily  global  solar  radiation  was  proposed  by 
Angstrom  [3],  who  derived  a  linear  relationship  between  the  ratio 
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Nomenclature 

So 

monthly  average  maximum  possible  daily  sunshine 

a-t 

regression  coefficients 

Sn/i 

duration  (h) 

sunshine  duration  that  take  into  account  the  natural 

c 

mean  total  cloud  cover  during  daytime  observations 

horizon  of  the  site  (h) 

(octa) 

ST 

mean  soil  temperature  (°C) 

CF 

cloud  factor 

T 

mean  air  temperatures  (  C) 

E 

mean  evaporation  (cm) 

Tk 

mean  air  temperature  (I<) 

H 

monthly  average  daily  global  radiation  on  horizontal 

T'max 

mean  maximum  daily  temperatures  (  C) 

surface  (MJ/m2  day) 

Train 

mean  minimum  daily  temperatures  (  C) 

Ho 

monthly  average  daily  extraterrestrial  radiation  on 

AT 

Temperature  difference  (°C) 

horizontal  surface  (MJ/m2  day) 

A  Tm 

monthly  mean  AT  (  C) 

Hc 

monthly  average  clear  sky  daily  global  radiation  on 

td 

mean  dew  point  temperature  (°C) 

horizontal  surface  (MJ/m2  day) 

V 

mean  water  vapor  pressure  (hpa) 

Hb 

monthly  average  daily  beam  radiation  on  horizontal 

W 

atmospheric  precipitable  water  vapor  per  unit 

Hi 

surface  (MJ/m2  day) 

monthly  average  daily  diffuse  radiation  on  horizontal 

z 

volume  of  air  (cm) 
altitude  of  site  (km) 

surface  (MJ/m2  day) 

es(T) 

saturation  vapor  pressure  at  the  air  temperature 

^mod 

the  estimated  global  solar  radiation  corrected  for 

T (kPa) 

f^day 

systematic  bias  (MJ/m2  day) 
half-day  length 

Creek  letters 

h 

k 

hourly  beam  radiation  (MJ/m2  h) 
hourly  diffuse  radiation  (MJ/m2  h) 

S 

solar  declination  (°) 

Isc 

solar  constant  (1367  W/m2) 

9 

solar  zenith  angle  (°) 

n 

number  of  day  of  year  starting  from  the  first  of 

V 

latitude  of  the  location  (°) 

January 

L 

latitude  of  the  location  (rad) 

p 

mean  precipitation 

X 

longitude  of  the  location  (°) 

Po 

atmospheric  pressure  at  sea  level  (hpa) 

cos 

sunset  hour  angle  (°) 

Ps 

atmospheric  pressure  at  the  site  (hpa) 

K 

the  noon  solar  altitude  angle  of  the  sun  (°) 

PS 

ratio  between  mean  sea  level  pressure  and  mean 

T 

atmospheric  transmittance 

daily  vapor  pressure 

^t.max 

maximum  (cloud-free)  daily  total  transmittance  at  a 

PWV 

precipitable  water  vapor  (mm) 

location 

RH 

mean  relative  humidity  (%) 

max 

proportion  of  Tt  max  observed  on  a  given  day  (cloud 

Rtimax 

maximum  relative  humidity  (%) 

correction) 

s 

monthly  average  daily  bright  sunshine  duration  (h) 

of  average  daily  global  radiation  to  the  corresponding  value  on  a 
completely  clear  day  (H/Hc)  at  a  given  location  and  the  ratio  of 
average  daily  sunshine  duration  to  the  maximum  possible  sun¬ 
shine  duration. 


H 

—  =a  +  b 


(1) 


A  basic  difficulty  with  Eq.  (1)  lies  in  the  definition  of  the  term 
Hc.  Prescott  [4]  and  the  others  have  modified  the  method  to  base 
it  on  extraterrestrial  radiation  on  a  horizontal  surface  rather  than 
on  clear  day  radiation  and  therefore  proposed  the  following 
relation: 


H 

Ho 


=  a+b 


(2) 


The  values  of  the  monthly  average  daily  extraterrestrial 
irradiation  (H0)  can  be  calculated  from  the  following  equation  [5]: 


71  ' 


/360n\ 

s(w) 


1  +  0.033  cos^ 

^coscp  cosS  sinms+  j^qCOs  sirup  s/nbj 


(3) 


The  solar  declination  (<5)  and  the  mean  sunrise  hour  angle  (ms) 
can  be  calculated  by  following  equations  [5]: 


S  =  23.45  sin 


'360 

365 


(284  +  n) 


(4) 


(5) 


For  a  given  month,  the  maximum  possible  sunshine  duration 
(monthly  average  day  length,  S0)  can  be  computed  by  using  the 
following  equation  [5]: 

S°=^ffls  (6) 

Although  the  Angstrom-Prescott  [4]  equation  can  be  improved  to 
produce  more  accurate  results,  it  is  used  as  such  for  many  applica¬ 
tions.  Some  of  the  regression  models  based  on  the  Angstrom-Prescott 
model  which  proposed  in  literature  are  given  as  follows: 

-  Jain  model  for  Italy  [6] 

^  =  0. 177  + 0.692  (+T)  (7) 

no 

-  El-Metwally  model  for  Egypt  [7] 

^-  =  0.228  +  0.527^  (8) 

no  V^o  / 

-  Bakirci  model  for  Turkey  [8] 

^  =  0.2786 +  0.41 60  (fj  (9) 

no  V-W 

-  Alsaad  model  for  Amman,  Jordan  [9] 

",0.174+0.6,5  (|) 


cos  =  cos  1  [—tamp  tand) 


(10) 
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Jain  and  Jain  model  for  Zambian  [10] 

^=  0.240  +  0.513  f/j  (11) 

no  VO/ 

Katiyar  et  al.  model  for  India  [1] 

=0.2281  +  0.5093  (f-j  (12) 

no  Vo  / 

Lewis  model  for  state  of  Tennessee,  U.S.A.  [11] 

|  ,0,4+0.57(|)  03, 

Tiris  et  al.  model  for  Turkey  [12] 

»,  0., 8  +  0.62(1)  (.4, 

Almorox  and  Hontoria  model  for  Spain  [13] 

^  =0.2170+0.5453  f/j  (15) 

no  VO  / 

Raja  and  Twidell  model  for  Pakistan  [14,15] 

=0.335  +  0.367  f  /  j  (16) 

no  V0  / 

Li  et  al.  model  for  Tibet,  China  [16] 

=0.2223  +  0.6529  ( /- ]  (17) 

no  \So  J 

Said  et  al.  model  for  Tripoli,  Libya  [17] 

4-  =  0.215  +  0.527(V)  (18) 

no  V0  / 

Ulgen  and  Ozbalta  model  for  Izmir-Bornova,  Turkey  [18] 

=0.2424  +  0.5014/^  (19) 

no  V0  / 

El-Sebaii  et  al.  model  For  Jeddah,  Saudi  Arabia  [19] 

"—2^, +3.78(|)  (20, 

Rensheng  et  al.  model  for  86  stations  in  China  [20] 
^=0.176  +  0.563^  (21) 

Ahmad  and  Ulfat  model  for  Karachi,  Pakistan  [21] 

4-  =0.324  +  0.405  (-p\  (22) 

no  \^0/ 

Jin  et  al.  model  for  69  stations  in  China  [22] 

=0.1332  +  0.6471  (f)  (23) 

Ho  \->o/ 


Akpabio  and  Etuk  model  for  Onne  region  (within  the  rainforest 
climatic  zone  of  southern  Nigeria)  [23] 

"  ,0.23+0.38(1)  ,24, 


Ulgen  and  Hepbasli  model  for  Ankara,  Istanbul  and  Izmir  in 


Turkey  [24] 

^  =  0.2671+  0.4754  (f)  (25) 

no  V>0/ 

Aras  et  al.  model  for  twelve  provinces  in  the  Central  Anatolia 
Region  of  Turkey  [25] 

=  0.3078 +  0.41 66  (26) 

Ho  \jo  J 

Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 

£-=  0.318  +  0.449  (f)  (27) 

no  V0/ 

Kholagi  et  al.  model  for  three  different  stations  in  Yemen  [27] 
^-  =  -0.191+0.571^0  (28a) 

^=  0.297  + 0.432  f/j  (28b) 

no  V0  / 

^  =0.262  + 0.454  (V)  (28c) 

Ho  v->o/ 

Katiyar  et  al.  model  for  Jodhpur,  Calcutta,  Bombay,  Pune  in 
India,  respectively  [1] 

^  =0.2276  +  0.5105  (+T)  (29a) 

no  VO/ 

=  0.2623  + 0.3952  (V)  (29b) 

no  VO  J 

^  =0.2229  +  0.5123  (V)  (29c) 

no  V0  J 

^  =  0.2286 +  0.5309  (f-|  (29d) 

Ho  \So  / 

Jain  model  for  (Salisbury,  Bulawayo  and  Macerata,  Italy), 
respectively  [28] 

=  0.3 13  +  0.474  (V)  (30a) 

no  V0  / 

=0.307  +  0.488  (f)  (30b) 

no  V^O/ 

0.309  + 0.599  Q0  (30c) 

Veeran  and  Kumar  model  for  two  tropical  locations  Madras 
and  Kodaikanal  in  India,  respectively  [29] 

^,0.34+0.32(|.)  (31a, 

^-=0.27 +  0.65  Q0  (31b) 

Chegaar  and  Chibani  model  for  Algiers  and  Oran,  Beni  Abbas 
and  Tamanrasset,  Algeria  [30] 

=  0.309  +  0.368  (f)  (32a) 

no  V^O  / 
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=  0.367  + 0.367 

no  V>0 

£-=  0.233  +  0.591  (j- 
no  V>0 


(32b) 

(32c) 


-  Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah  weather 
stations  in  Oman,  respectively  [31] 


=  0.3326  +  0.3110("  rS 

no  V>0 

^  =0.2418  +  0.3555^ 

no  V>0 


(33a) 

(33b) 


2.1.2.  Model  2:  Glower  and  McCulloch  model 

Glower  and  McCulloch  proposed  the  following  equation, 
which  takes  into  account  the  effect  of  latitude  of  the  site,  cp,  as 
an  additional  input  and  is  valid  for  ( cp  <  60°) [32]: 


A,  ,(! 


^  =0.29  coscp  +  0.52  ( ^ 

no  V>0 


(34a) 


Some  other  corresponding  modified  relations  are: 


-  Ulgen  and  Hepbasli  model  for  Izmir,  Turkey  [33] 

=  0.3092  cosy) +  0.4931  ( £ 
no  V^O 


-  Raja  and  Twidell  model  for  Pakistan  [14,15] 
0  =  0.388  coscp  +  0.367  ( ^ 

no  V>0 


^  =  0.388  coscp  +  0.407  ( f 
no  V^O 


(35) 

(36a) 

(36b) 


2.1.3.  Model  3:  Page  model 

Page  has  provided  the  coefficients  of  the  modified  Angstrom- 
type  model  [4],  which  is  claimed  to  be  applicable  anywhere  in  the 
world  [34]: 


wr023+0AH 


(37) 


2.1.4.  Model  4:  Rietveld  model 

Rietveld  by  using  measured  data  collected  from  42  stations 
located  in  different  countries,  has  proposed  a  unified  correlation 
to  compute  the  horizontal  global  solar  radiation.  Rietveld’s  model, 
which  is  claimed  to  be  applicable  anywhere  in  the  world,  is  given 
in  following  equation  [35]: 


tJ-  =  0.18  +  0.62  Tir- 

no  V^O 


(38) 


Rietveld  also  examined  several  published  values  of  a  and  b 
coefficients  of  the  Angstrom-Prescott  model  [4]  and  noted  that 
constants  a  and  b  are  related  linearly  to  the  appropriate  mean 
value  of  S/So  as  follows  [35]: 

'S' 


a  =  0.10  +  0.24 


b  =  0.38  +  0.08 


So 


(39a) 


(39b) 


2.1.5.  Model  5:  Dogniaux  and  Lemoine  model 

Dogniaux  and  Lemoine  proposed  the  following  correlation, 
where  the  regression  coefficients  of  the  Angstrom-Prescott  model 
[4]  seem  to  be  as  a  function  of  the  latitude  of  the  site  [36]: 


H 

Ho 


=  0.37022  + 


0.005061  j-  [-0.00313 


cp  +  0.32029 


(40) 


(34b)  April 


They  also  in  same  year  obtained  the  specific  monthly  correla¬ 
tions  which  are  listed  in  Eq  (41a)-(411)  [36]: 

January  //  =  (-0.00301c/>  +  0.34507)  +  (0.00495y)  + 0.34572) 

n  0  VO/ 

(41a) 

February  =  (-0.00255c/)  +  0.33459)  +  (0.00457 cp  +  0.35533)  ( ^ 

no  \S(), 

(41b) 

March  =  (-0.00303c/)  +  0.36690)  +  (0.00466</)  + 0.36377)  [^- 

no  V>0, 

(41c) 
S' 


=  (-0.00334c/)  +  0.38557)  +  (0.00456c/)  +  0.35802) .  ^ 


May 


June 


H 

Ho 


=  (-0.00245c/) +  0.35057)  +  (0.00485c/)  + 0.33550) 


(41  d) 
S' 


VSo. 

(41  e) 


//  =  (-0.00327</)  + 0.39890)  +  (0.00578y)  +  0.27292)  ( f 

no  V>0 


(41  f) 


July  —  =  (-0.00369y)  + 0.41 234)  +  (0.00568y)  + 0.27004)  ( — 
Ho  \5o 


(41  g) 


August  F-  =  (-0.00269c/) +  0.36243)  +  (0.00412y)  +  0.33162) 

no  V^0 


(41  h) 


H 


September  =  (-0.00338c/)+0.39467)  +  (0.00564c/)  +  0.27125) 


H 


(41  i) 
S 


October  =  (-0.00317</)  +  0.36213)  +  (0.00504c/)  +  0.31790) . 

H  o  \So 


November  "  =  (-0.00350</)  + 0.36680)  +  (0.00523y)  +  0.31467) 


(41  j) 
S 


(41k) 

S' 


December  —  =  (-0.00350</>  + 0.36262)  +  (0.00559c/>  +  0.30675)  ( - 
no  \So. 


(411) 


2.1.6.  Model  6:  Kilic  and  Ozturk  model 

Kilic  and  Ozturk  calculated  the  a  and  b  regression  coefficients 
of  Angstrom  model  [4]  for  Turkey  [37]: 


a  =  0.103  +  0.000017  Z+0.198cos(c/>-<5) 
b  =  0.533-0.165  cos(cp-d) 


(42a) 

(42b) 
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2 A. 7.  Model  7:  Benson  et  al.  model 

Benson  et  al.  proposed  two  different  formulations  for  two 
intervals  of  a  year  depending  on  the  climatic  parameters  to 
estimate  the  global  solar  radiation  [38]: 

For  January-March  and  October-December 


|-0,8+0<l 


For  April-September 


A. 0.24  +  0.53(1 


(43a) 


(43b) 


2.3.8.  Model  8:  Ogelman  et  al.  model 

Ogelman  et  al.  expressed  the  ratio  of  global  to  extraterrestrial 
radiation  by  a  second  order  polynomial  function  of  the  ratio  of 
sunshine  duration  [39]: 


*-‘©«(s 


m  ■0  '95  +  0-676(I)— °  '42(| 


(44a) 


(44b) 


Many  authors  in  all  over  the  world  applied  this  model  and 
determined  the  regression  coefficient  of  this  model  for  particular 
location  of  interest  as  follows: 


-  Akinoglu  and  Ecevit  model  for  Turkey  [40] 


=0.145  +  0.845fJ-V°-280f  c 

no  v>o  /  V^o. 


-  Almorox  and  Hontoria  model  for  Spain  [13] 


=0.1840  +  0.6792 
no 


(I  )-01228(: 


!)' 


-  Bakirci  model  for  Turkey  [8] 


=0.2545  +  0.5121  (©-] -0.0864]  ^ 
no  \Sq  J  W 


(45) 


(46) 


(47) 


-  Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 

0.1541  +  1.17, 4(|)-0.705(*)'  (53) 

-  Ahmad  and  Ulfat  model  for  Karachi,  Pakistan  [21] 

.  2 


„  =  0.348  +  0.320^  CS^1  +  0.07of  ^ 
no  \^0/  V^o 


(54) 


-  Tahran  and  Sari  model  for  Central  Black  Sea  Region  of  Turkey  [44] 


H 

Ho 


=  0.1874  +  0.8592 


©-0-4764©1 


-  Jin  et  al.  model  for  69  stations  in  China  [22] 


^  =  0.1404  +  0.6126(0  +0.0351  (^- 


(55) 


(56) 


-  Aras  et  al.  model  for  twelve  provinces  in  the  Central  Anatolia 
Region  of  Turkey  [25] 

^  =  0.3398  +  0.2868  (JPj  +0.1187  (jPj  (57) 

-  Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah  weather 


stations  in  Oman,  respectively  [31] 

=0.9428-1.2027  ( +0.9336^  Ap) 
no  V^0  /  vW 

=0.1971  +0.6297  (£-) -0.2637  fp-) 
Ho  v->o/  \jo / 


2.3.9.  Model  9:  Bahel  et  al.  model 

Bahel  et  al.  reported  the  following  relationship  [45]: 

^  =  0.1 75  + 0.552 
Ho  \5q 


(58a) 

(58b) 


(59) 


Tasdemiroglu  and  Sever  model  for  Ankara,  Antalya,  Diyarbakir, 
Gebze,  Izmir  and  Samsun  in  Turkey  [41] 

«  ,0.225  +  0.014(1)  +0.001  ,48, 

Yildiz  and  Oz  model  for  Turkey  [42] 

^  =0.2038 +  0.9236  f^V0-391  (49) 

Ho  \5o/  \5o  / 

Aksoy  model  for  Ankara,  Antalya,  Samsun,  Konya,  Urfa  and 
Izmir,  Turkey  [43] 

^=0.148  +  0.6680-0.0790  (50) 

Said  et  al.  model  for  Tripoli,  Libya  [17] 

£-0-,+a874©-“-255©! 

Ulgen  and  Ozbalta  model  for  Izmir-Bornova,  Turkey  [18] 

=  0.0959  +  0.9958  -0.3922  f (52) 

Ho  \jo  /  V>oJ 


2 A A0.  Model  30:  Zabara  model 

Zabara  correlated  monthly  a  and  b  values  of  the  Angstrom- 
Prescott  model  [4]  with  monthly  relative  sunshine  duration  (S/S0) 
as  a  third  order  function  and  expressed  the  a  and  b  coefficients  as 
[46]: 

a  =  0.395-1.247  0  +2.680  (JPj  -1.674 0  (60a) 

b  =  0 .395  +  1. 384  (0j  -3.249  ^0  +  2.0550  (60b) 

2.3.3  3.  Model  3  3;  Bahel  model 

Bahel  developed  a  worldwide  correlation  based  on  bright 
sunshine  hours  and  global  radiation  data  of  48  stations  around 
the  world,  with  varied  meteorological  conditions  and  a  wide 
distribution  of  geographic  locations  [47]. 

A  =  0.16+0.87  (7)  —0.6,  (7)  2 +0.34(1-)  3  <61, 

2.3.32.  Model  12:  Copinathan  model 

Copinathan  suggested  a  and  b  regression  coefficients  of 
Angstrom-Prescott  model  [4]  as  a  function  of  elevation  (Z)  and 


F.  Besharat  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  2?  (2013)  798-821 


805 


sunshine  ratio  (S/S0)  for  estimation  of  the  global  solar  radiation  [48]: 

S  '  (62a) 


a  =  0.265  +  0.07Z— 0.135 


b  =  0.401— 0.108Z+0.325 


(62b) 


H 

Ho 


Gopinathan  also  reported  the  following  correlations  [49]: 

S' 


-0.309  +  0.539  cos<p— 0.0693Z + 0.290 


+ 


1.527-1.027  coscp  +  0.0926Z- 


°-359(I) 


(63) 


LJ 

September  —  =  0.20+0.59 
no 

October  ©  =  0.19  +  0.60  (© 


LJ 

November  —  =  0.17  +  0.66 

no 


LJ 

December  —=0.18  +  0.65 
Ho 


(671) 

(67j) 

(67k) 

(671) 


Also,  they  reported  the  regression  coefficients  of  a  and  b  as 
follows  [51]: 


2.1  A3.  Model  13:  Newland  model 

A  linear-logarithmic  model,  has  been  used  by  Newland  to 
obtain  the  best  correlation  between  H/H0 and  S/S0  [50]: 


H 

Ho 


H 


=  a+b 


+  clog 


Ho 


—  =0.34  +  0.40  —  +0.171og(  — 


So 


So, 


(64a) 

(64b) 


Some  applications  of  this  model  are: 

-  Bakirci  model  for  Turkey  [8] 

©  =  0.3925  +  0.2877)©)  +0.1527  logl  f 


Ho 


So 


So 


(65) 


-  Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah  weather 
stations  in  Oman,  respectively  [31] 


©  =  -1.1931  +  1.8641  (  f- )  -1.2544  log[  f- 


©  =0.4337  +  0.1430 
no 


1) 


+  0.0861  log  [  — 


(66a) 

(66b) 


a  =  0.179  +  0.099 U- 


b  =  0.1640  +  0.1786 


-,-0935(l) 


(68a) 

(68b) 


2.1.15.  Model  15:  Luhanga  and  Andringa  model 

Luhanga  and  Andringa  derived  their  own  model  as  follow  [52]: 


=0.241+0.488 

no 


(I; 


(69) 


2.1.16.  Model  16:  Louche  et  al.  model 

Louche  et  al.  presented  the  model  below  to  predict  global  solar 
radiation  [53]: 


H 

Ho 


=  0.206  +  0.546 


(70) 


Furthermore,  the  Angstrom-Prescott  model  [4]  has  been 
modified  through  the  use  of  the  ratio  of  (S/Snh)instead  of  (S/S0) 
by  Louche  et  al.  model  which  is  presented  as  follows  [53]: 


2.1.14.  Model  14:  Monthly  specific  Rietveld  model  (Soler  model) 
Soler  applied  Rietveld’s  model  to  100  European  stations  and 
gave  the  following  specific  monthly  correlations  [51]: 


LJ 

January  —  =  0.18  +  0.66 
no 

LJ 

February  —=  0.20+0.60 
no 


March  ^=0.22  +  0.58 
no 


April  ^=0.20  +  0.62^ 


May 


Ho 

H 

Ho 

H 


=  0.24  +  0.52 


June  —  =0.24+ 0.53 


July 


Ho 

H 

Tfa 


=  0.23  +  0.53 


H 


August  —=0.22  +  0.55)^ 


(67a) 

(67b) 

(67c) 

(67d) 

(67e) 

(67f) 

(67g) 

(67h) 


H 

Hi 


=  a+b 


o 

1 

Sn(i 


s 

Snh 


0.8706 


->o 


+  0.0003 


(71a) 

(71b) 


2.1.17.  Model  17:  Samuel  model 

Samuel  has  correlated  H/H0with  S/S0  in  the  form  of  a  third 
order  polynomial  equation  [54]: 


H 

Ho 

H 

Ho 


=  a+b 


+c 


=  -0.14  +  2.52 


+  d 


-3.71 


So 


+  2.24 


(I)’ 


(72a) 


(72b) 


Some  investigator  used  this  model  to  estimate  global  solar 
radiation  in  different  locations  all  over  the  world.  Some  of  these 
models  are: 


-  Ertekin  and  yaldiz  model  for  Antalya,  Turkey  [55] 


=-2.4275  + 11.946)  ^  ) -16.745)  ^  )  +7.9575)  ^ 


© 


© 


(73) 
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-  Almorox  and  Hontoria  model  for  Spain  [13] 

0=0.230  +  0.3809000.46940)  -0.3657  0)  (74) 

-  Lewis  model  for  state  of  Tennessee,  U.S.A.  [11] 

wrom-33il)+73iff-4  Hff  <751 

-  Ulgen  and  Hepbasli  model  for  Izmir,  Turkey  [33] 

0  =  0.2408+0.3625000.45970)  -0.37080)  (76) 

-  Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 

0=0.1796+  0.9813  000.2958  0)  -0.2657(0  (77) 

-  Ulgen  and  Hepbasli  model  for  Ankara,  Istanbul  and  Izmir  in 
Turkey  [24] 

0  =  0.2854+0.25910  +0.6171  0)  -0.48340  (78) 

-  Tahran  and  Sari  model  for  Central  Black  Sea  Region  of 
Turkey  [44] 

"=0.1520+1.13340— 1.112602+  O.451603  (79, 

-  Jin  et  al.  model  for  69  stations  in  China  [22] 

0=0.1275+0.7251000.22990)  +0.18370  (80) 

-  Aras  et  al.  model  for  twelve  provinces  in  the  Central 
Anatolia  Region  of  Turkey  [25] 

0  =0.4832-0.6161  (0)  +1.8932(0)  -1.0975  ("0 
"0  \J0  /  V>0/  VO/ 

(81) 

-  Rensheng  et  al.  model  for  86  stations  in  China  [20] 

A  =  0.150+1.145(1-)— 1. 474(A)2 +o.963(T)  3  (82) 

-  Bakirci  Model  for  Erzurum,  Turkey  [56] 

0  =  0.6307-0.7251  001. 20890 )  -0.4633  0)  (83) 


2.3.18.  Model  18:  Gopinathan  and  Soler  model 

Gopinathan  and  Soler  suggested  following  linear  equations  for 
locations  with  latitudes  between  60°N  and  70  N  [57]: 

0=0.1538  +  0.7874(0)  (84a) 

no  V^O  / 


0  =  0.1961+0.7212f  0 
no  VW 


(84b) 


2.3.39.  Model  39:  Raja  model 

Based  on  Bennett's  formula  [58],  Raja  proposed  the  following 
insolation-sunshine  relation  [59]: 


H 

Ho 


0.368-0.125  1- 


0.667-0.01 8  (  1  — —  )  —0.21 1  coscp 


S 

So4 


(85a) 


where  Za  =  8000  in  and  SM  is  the  4°  corrected  day  length  used  to 
compensate  the  finite  threshold  of  the  Champbell-Stokes  sunshine 
recorder  and  is  calculated  as  follows: 


S04  = 


2  _j  r(sin4°-sin<psin<5) 
15  [  cosdcoscp 


(85b) 


2.3.20.  Model  20:  Coppolino  model 

Coppolino  developed  a  power  function  and  incorporated  a 
trigonometric  term  to  estimate  global  solar  radiation  [60]: 


(a)  Power  model 


(86) 


-  Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah 
weather  stations  in  Oman,  respectively  [31] 


H 

Ho 


£-0.4470 


(87a) 


H 

Ho 


(b)  Power-trigonometric  model 


H 

Ho 


So 


( sinhaf 


(87b) 


(88) 


The  noon  solar  altitude  angle  of  the  Sun  (ha)  is  the  comple¬ 
ment  of  zenith  angle  and  can  be  calculated  from  [60]: 

hx  =  90-cp+d  (89) 

-  Coppolino  employed  this  model  to  estimate  global  solar 
radiation  in  34  Italian  stations.  The  regression  coefficients 
of  this  model  are  presented  as  follows  [60] 

0=0.670)  ( sinK )0  05  0.15  00  <0.90  (90) 

-  Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah 
weather  stations  in  Oman,  respectively  [31] 

i_i  /  C  \  0.5239 

0=e-O41350)  (sinhx)009n  (91a) 

0=e-O  559S0)  (sinha)-00296  (91b) 

2.3.23.  Model  23:  Tin's  et  al.  model 

Tiris  et  al.  reported  their  own  correlations  as  follow  [61]: 

0  =0.2262  +  0.41 8  (0)  (92) 

Ho  \SoJ 

2.3.22.  Model  22:  Ampratwum  and  Dorvlo  model 

Ampratwum  and  Dorvlo  have  suggested  a  logarithmic  type 
model  as  [31]: 

wra+blog(£) 


(93) 
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Some  other  developed  logarithmic  correlations  are: 

-  Almorox  and  Hontoria  model  for  Spain  [13] 

=0.6902  +  0.6142  log( f 
no  \V 


(94) 


Ampratwum  and  Dorvlo  model  for  Seeb  and  Sallalah  weather 
stations  in  Oman,  respectively  [31] 


D  =  0.6376  +  0.2490  log( ^ 
no  \$0. 

=0.5612  +  0.1412  log(^ 
no  V^O 


(95a) 

(95b) 


a  =  — 0.0344/n  (j)  +0.1982;  b  =  -0.0201  In 


a  =  0.1950  +  0.2943  (j-j  -0.3706  (J- 


(!) 


+  0.4562 
(101b) 


b  =  0.6248-0.7033 


I  +0.6368 


a  =  0.2646- 


b  =  0.5403  +  0.0149 


(!) 

°-2958(!)+0-967,(!)!-°-8749(!) 
;)-°",3(!)2+,  0649(!) 


(101c) 


(lOld) 


-  Bakirci  model  for  Turkey  [8] 

=0.6446  +  0.4842  log( $ 
no  \V 


(96) 


Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 


^  =0.698  +  0.2022  log(pj 


(97) 


a  =  0.2226  +  0.2841  (pj  -1.3653  (p 

+2.6661  (|)Y, 7"7(| 


b  =  0.561—0.2715 


(!) 


+  0.1606 


So 


-°-6839(!) 


—  )  +0.8888 


(lOle) 


2.1.23.  Model  23:  Klabzuba  et  al.  model 

Klabzuba  et  al.  developed  the  following  relationship  between 
the  measured  H  and  the  relative  sunshine  duration  [62]: 


H  =  a+b  |+c 


d+ 


(n-e)2 


(98) 


-  Klabzuba  et  al.  model  for  Hradec  Kralove  solar  observatory, 
Czech  Republic  [62] 

H=  7.19  +  0.258^^—9.28  x  10~6 

2.1.24.  Model  24:  Togrul  et  al.  model 

Togrul  et  al.  applied  the  various  regression  analyses  to  inves¬ 
tigate  the  relation  between  monthly  mean  S/S0  ratio  and  con¬ 
stants  a  and  b  in  Angstrom  model  [4]  as  follow  [63]: 

a=i+j'(£);  b=k+i(l)  <iooa) 

a  =  iln(pj+j ;  b  =  kln(pj+m  (100b) 


So  +22'9 


(n-174.7)2  (99) 


-  Togrul  et  al.  obtained  the  empirical  coefficients  of  their 
models  for  Elazig,  Turkey  as  follows  [63] 

a  =  0.281 6-0.1 049  (jpj :  b  =  0.4762-0.01 74  (pj 

(101a) 


2.1.25.  Model  25:  Elagib  and  Mansell  model 

Elagib  and  Mansell  developed  new  techniques  for  predicting 
solar  radiation  based  on  sunshine  hours  and  geographical  para¬ 
meters  [64]: 


H 

Ho 


=  a  exp 


(102a) 


^-=a+bL+cZ+d(f)  (102b) 

Ho  V-W 


H 

Ho 


a+bL+c 


(102c) 


=  a+bZ+c(-^)  (102d) 

Ho  \So  / 

Some  researchers  calibrated  the  Elagib  and  Mansell  models  for 
different  locations.  Some  of  these  models  are  as  follows: 


-  Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 


H 

Ho 


0.3396  exp 


0.8985 


(103) 


-  Jin  et  al.  model  for  69  stations  in  China  [22] 

(I  =  0. 0855  +  0. 0021 +  0.03Z+0. 5654  f  (S')  (104) 

no  V0/ 


-  Rensheng  et  al.  model  for  86  stations  in  China  [20] 

=0.122  +  0.00U  +  2.57  x  10~2Z  +  0.543  (£-)  (105) 

Ho  \5o  J 


2.1.26.  Model  26:  Monthly  specific  Elagib  and  Mansell  model 

Elagib  and  Mansell  have  investigated  the  possibility  of  estab¬ 
lishing  monthly-specific  equations  for  estimating  global  solar 
radiation  across  Sudan.  The  best  performing  equations  for  each 
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month  are  given  as  following  [64]: 

January  =  0.1357  +  03204L  +  0.0422Z+0.4947  ( 

no  \-->o 

February  =  0.1563  +  0.3166L  +  0.1006Z+0.4593(" 

n  o  V>0 


March  =  0.7727 

no  V^O 


April  =  0.1640  +  0.0397Z+0.5773(^ 

no  \V 


May  =0.0709  +  0.8967  (Jpj -0.2258  (J- 


June  ^  =  -0.0348  +  1.5078^-0.8246^ 


July  =0.3205  +  0.1444L  +  0.0782Z+0.2916(/^ 

no  W 


August  =  0.2720  +  0.0369L+0.1017Z+0.3888( 

n  o  v>0 


September  ^=-0.3710+2.5783^-1.6788^ 


H 


il) 


October  —  =  0.1593-0.1043L  +  0.0609Z  +  0.5916(  —  ] 


November  =  0.1786  +  0.0199Z+0.5441 
no 


December  =  0.1714  +  0.1329L  +  0.0482Z  +  0.5015(/  ^ 

no  \J0 


(106a) 

(106b) 

(106c) 

(106d) 

(106e) 

(106f) 

(106g) 

(106h) 

(106i) 
( 1 06j ) 
(106k) 

(1061) 


2.2.27.  Model  27:  Togrul  et  al.  model 

Togrul  et  al.  developed  some  statistical  relations  to  estimate 
monthly  mean  daily  global  solar  radiation  by  using  clear  sky 
radiation.  The  correlations  have  been  developed  by  employing 
both  ratios  of  (S/S0)  and  (S/Snh)[65J. 

where  Snh  is  the  monthly  mean  sunshine  duration  that  take 
into  account  the  natural  horizon  of  the  site  and  can  be  calculated 
by  following  equation: 

1  =0^06+00003  (!07) 
5nJi  50 

For  summer 
H 


Hc 


=  0.5771+0.6843 


(!)-» 


-0.3544 


=  0.4356  + 1.3939  (^[-1.4988 
nc  \JQ 


(l): 


+  0.5954 


=  0.7295exp 


0.2633 


H 
He 

=0.9337.  c 
nc  W0 


=0.577  +  0.7825 (  ^-Vo.4631(/- 

n  c  \^nhj  \^nh 


=0.4356+ 1.5935  ( -1.9584  f- +0.8893^ 
nc  \JnhJ  \JnhJ  \^nh 


£-°-9552(£> 

For  winter 


=0.2384  +  2.1498 ( (T)-1.8664(  $ 
nc  W0/  WO. 


^.03081  +  1.4933(|)-0.0839(|-)2-,. 4513(|)3 


=  0.5451  exp 

(I) 


0.7981 


=  1.0302 


H 
Hc 

H 

He 

£-=  0.2385  +  2.4604  fj-') -2.445  ( J~ 

nC  Wnh/  \^nh 


H 


=  0.3077  +  1.7137(^-0.1251  ( r~)  -2.1612  (■ 
nc  \JnhJ  \Jnh/  \^nhj 


(108f) 

(108g) 

(109a) 

(109b) 

(109c) 

(109d) 

(109e) 

(109f) 

(109g) 


2.1.28.  Model  28:  Almorox  and  Hontoria  model 

Almorox  and  Hontoria  proposed  the  following  exponential  corre¬ 
lation  to  estimate  global  solar  radiation  from  sunshine  hours  [13]: 


l-+bexp(i 


-  Almorox  and  Hontoria  model  for  Spain  [13] 


=  -0.0271  +0.3096  exp ( £ 

no  W0. 


-  Bakirci  model  for  Turkey  [8] 

=0.0963  +  0.2337  exp  ( ^ 
no  W0 


(110) 


(ill) 


(112) 


2.1.29.  Model  29:  Jin  et  al.  model 

Jin  et  al.  proposed  the  following  models  by  using  solar 
radiation  data  and  some  geographical  parameters  like  latitude 
and  altitude  [22]: 


(108a) 

=a+bcos<p+cZ+d(Jpj 

(113a) 

(108b) 

d  =  (a+bL+cZ)+(d+eL+fZ)l 
H  0  ' 

I) 

(113b) 

(108c) 

—  =  (a+bcoscp+cZ)  +  (d+ecoscp+fZ)  | 

f£) 

(113c) 

(108d) 

^  =(a+bL+cZ)  +  (d+eL+fZ)l 
no  ' 

+)+(g+M+a(|.)2 

(113d) 

(108e) 

—  =  (a+bcos(p+cZ)  +  (d+ecoscp+fZ)  | 

[Q 

+  (g+hcoscp+iZ)  (JPj 
Some  implications  of  these  models  are: 

-  Jin  et  al.  model  for  69  stations  in  China  [22] 
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H 


(113e)  April 
May 


H 

Ho 

H 

Ho 


=  2.1 186-2.0014  cos(p  +  0.0304Z  +  0.5622(  — 


=  (0. 1094  +  0. 0014L+0.0212Z) 

+  (0.5176  +  0.0012L+0.0150Z) 


=  (1.8790-1.7516  cos(p  +  0.0205Z) 
Ho 

+  (1.0819-0.5409  coscp  +  0.01 69Z) 
=  (0.021 8  +  0.00331  +  0.0443Z) 


(114a) 


(114b) 


(114c) 


H 

Ho 


+  (0.9979— 0.00921— 0.0852Z) 


H 

Ho 


+  (— 0.5579 +  0.01 20(p  + 0.1 005Z)  (^~ 


=  (4.2510-4.1878  C0S(?  + 0.0437Z) 

+  (-10.5774+1 1 .4512  coscp-0.0832Z) 


(114d) 


(114e) 


-  Rensheng  et  al.  model  for  86  stations  in  China  [20] 

=0.280-0.141  costp  +  2.60  x  10~2Z  +  0.542 (jPj  (115a) 


Ho 

H 

Ho 


=  ^0.329—0.196  coscp  +  2.2  x  10~2ZJ 
+  (^0.457  +  0.097  cos<p  +  6.72  x  10 ~3Z 


(115b) 


2.1.30.  Model  30:  El-Metwally  model 

El-Metwally  proposed  a  non-linear  correlation  between  clear 
index  (H/H0)  and  relative  sunshine  (S/S0)  as  follows  [7[: 

H 


H0 


=  a(ys/s0) 


-  El-Metwally  model  for  Egypt  [7] 
H 


—  =  0.7 1 3(l/s/s°) 
Ho 


(117) 


2.1.31.  Model  31:  Almorox  et  al.  model 

Almorox  et  al.  reported  the  monthly-specific  equations  for 
estimating  global  solar  radiation  from  sunshine  hours  for  Toledo, 
Spain  as  given  by  following  correlations  [66]: 

January  ^-=0.285  +  0.444 

H0  \So 


February  ^=0.272  +  0.465 
no 

March  ^- =  0.291 +0.491  (  ^ 


June 

July 


H0 

H 

Ho 

H 

Ho 

H 

Ho 


=  0.266  +  0.495 

=  0.286  +  0.475 

=  0.311+0.439 

=  0.329  +  0.406 
H 


S 

So 

s 

So 

s 

So 

s 

So 


August  —  =  0.313  +  0.410 
Ho 


September  ^  =  0.271  +0.479 
H0 

October  ^-=0.259  +  0.4651^ 


November  =  0.279+0.431  (  JJ- 
December  ^-=0.282  +  0.4281^ 


809 

(11 8d) 

(11 8e) 
(1180 
(H8g) 
(1 1 8h) 
(1 1  Si) 
( 1 1 8 j ) 
(118k) 
(1181) 


2.1.32.  Model  32:  Rensheng  et  al.  model 

Rensheng  et  al.  have  suggested  new  equations,  based  on  the 
Angstrom  model  [4]  to  estimate  global  radiation.  These  equations 
are  as  follows  [20]: 


H 

Ho 

H 

Ho 

H 

Ho 


H 

Ho 


.w+tL+cZ>+d(!) +e(!)M! 

.(a+to^+cZl+^D+eg)  +f 

=  (a+bcoscp  +  cZ)  +  (d+ecoscp+fZ)  (Jp\ 

+  (g + hcoscp  +  i Z)  (^JpJ  +  (/+  kcostp  +  IZ)  (Jpj 


—  (fl  +  bl+ cZ+  dl)  +  e  ( 


H 

Ho 


(116)  —  =  (a+bcos(p+cZ+dX)+e(J^j 

=  (a+bL+cZ+dA)+(e+fL+gZ+hA)t 


H 

Ho 


S 

VSo. 


H 

Ho 

H 

Ho 


=  (a + bcostp  +  cZ+  dl)  +  (e  +fcoscp+gZ + hi) 
=  (a+bL+cZ+dA+eA2)  +/ 


+h 


(118a,  ^  =  (a  +  bcos^  +  cZ+dl+el2)+/g)+g(£)  +KI 


(119a) 

(119b) 

(119c) 

(1 1 9d) 

(119e) 

(1190 
(1 19g) 

(1 19h) 

3 

(11 9i) 


(!) 

(118b) 

-  =  ( 

H0  V 

(118c)  +  ( f+gcoscp  +  hZ+iA+jA 2)  (Jpj 
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+  (k + Icoscp  +  mZ+nZ+ 0/ 12)  (pj 
+  (p+qcoscp+rZ+sl+tZ2J  (pj 


(H9j) 


Rensheng  et  al.  employed  their  models  to  obtain  following 
correlations  using  data  from  1994  to  1998  at  86  stations  in 
China  [20]. 

H 

TT0 


=  (0.109  +  0.001L  +  2.41  x  10~2z)  + 1.029  JJ 


H 

Ho 


-,.2.6(|)'+°.787(| 


=  (0.234—0.112  costp  +  2. 43  x  10~2Z 


+  1.026(|)-,. 209(|)  +0.782(1 


(120a) 


(120b) 


H 

Ho 


=  (0.336-0.233  coscp+2. 64  x  10-2z) 
+  (-0.670  +  2.140  coscp-OAZ)  (Jpj 

+  (2.744-5  coscp +  0.3Z)  (Jpj 


+  (-1.638  +  3.042  cos(p-0.2Z) 


H 

Ho 


H 

Ho 


=  (0.1 17— 0.001  L  + 2.59  x  10-2Z+4.11  x  lO^ZJ 
S' 


+  0.543 


So 


=  (0.275-0.141  coscp  +  2.63  x  10~2Z 


+  4.27  x  10  SA)  +0.542 


(120c) 


(120d) 


(120e) 


H 

Ho 


=  0.094  +  0.002L+2.27  x  10  2Z+ 0.0001  X 


+  ( 0.586— 0.0008L+ 5.36  x  1Q-3Z-0.00022 ) 


So 
(120f) 


H 


=  (0.313-0.195  cos<p  +  2.28  x  10~2Z+0.00012J 
+  (0.476  +  0.097  coscp+5.69  x  10~3Z-0.00022)  (pj 

(120g) 

=  (0.370  +  0.0007L  +  2.44  x  10~2Z-0.0052  +  2.24  x  10~522) 


— 0.1Z+ 0.0692— 0.003 22)  (pj  +  (l2.402-3.867  coscp 


+  0.3Z-0. 1992  +  0.00092" 


So 


-9.442 


+  2.115  cos(p-0.2Z+0.1642-0.0008Z2)  (p 


(120j) 


2.1.33.  Model  33:  Sen  model 

Sen  proposed  a  nonlinear  model  for  the  estimation  of  global 
solar  radiation  from  available  sunshine  duration  data.  This  model 
is  an  Angstrom  type  model  with  a  third  parameter  appears  as  the 
power  of  the  sunshine  duration  ratio  [67]: 


wr°+H 


(121) 


-  El-Sebaii  et  al.  calibrated  this  model  For  Jeddah,  Saudi  Arabia  [19] 

H  /  r  x  2.344 

^-  =  -0.864  +  1.862  -M  (122) 

Ho  \3o/ 


2.1.34.  Model  34:  Bakirci  model 

Bakirci  reported  the  original  Angstrom-type  equations  including 
the  linear,  second-order  and  fifth-order  polynomial  relationships 
between  the  monthly  average  values  of  (H/Hc)  and  (S/S0)  as 
follows  [68]: 


=0.7836-0.0460  ( 

He  V^o. 


£,..0.92— ..0547(1)  +  0.966,  (| 


(123a) 

(123b) 


P  =  -U. 225  + 128.010  (p -516.900  (pj  +  994.730^0 


-920.350  (p)  +329.93  (p 
v>o /  V^o 


(123c) 


2.1.35.  Model  35:  Bakirci  model 

Bakirci  developed  following  models  which  is  a  derivation  of 
the  modified  Angstrom-type  equation  [8[: 


(a)  Linear  exponential 

a  ,<,+„(!-) +«*,(! 


(124) 


1  026lsoJ-1-208UJ  +a783VV 


(120h) 


-  =  (o. 426— 0.087  coscp +2A4  x  10~2Z-0.0042  +  1.86  x  10“522) 

O'  ' 


(l)J+0-779(l 


H 

Ho 


- 1. 024  (^  1-1.204 


=  (1.012-0.141  cos<p  + 2.84  x  10-2Z-0.0142 


(120i) 


+  6.7  x  10~522)  +  (-4.061  +1.740  coscp 


Some  other  relations  are: 

-  Muzathik  et  al.  model  for  state  of  Terengganu,  Malaysia  [69] 

(125) 


H  =  0.19490  +  0.4771  (p)  +0.02994 exp  (  S 


H0 


So 


-  Bakirci  model  for  Turkey  [8] 


=0.3448  +  0.5636  ( p) -0.0838exp(  ) 
Ho  V^o /  V-W 


(126) 
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(b) Exponent 


-  Bakirci  model  for  Turkey  [8] 


H 

Ho 


=  0.6660 


(127) 


(128) 


-  Togrul  and  Togrul  model  for  Ankara,  Antalya,  Izmir,  Yenihisar 
(Aydin),  Yumurtalik  (Adana)  and  Elazig  in  Turkey  [26] 


H 

Ho 


=  0.7316 


(129) 


2.2.  Cloud-based  models 

Clouds  and  their  accompanying  weather  patterns  are  among  the 
most  important  atmospheric  phenomena  restricting  the  availability 
of  solar  radiation  at  the  earth’s  surface.  The  cloudiness  as  a  limiting 
factor  causes  distribution  and  dissipation  of  the  solar  radiation 
reaching  the  atmosphere  and  affects  the  amount  of  radiation 
received  at  the  earth’s  surface.  The  cloud  data  are  detected 
routinely  by  meteorological  satellites,  so  a  number  of  models  have 
been  developed  to  estimate  global  solar  radiation  from  observations 
of  various  cloud  layer  amounts  and  cloud  types  [70].  In  the 
following  section,  correlations  which  use  only  the  cloudiness  are 
presented  and  classified  based  on  their  developing  year. 

2.2.1.  Model  2:  Black  model 

Black  developed  following  quadratic  equation,  using  data  from 
many  parts  of  the  world  [71]: 

=0.803-0.340C-0.458C2  (C<0.8)  (130) 

Ho 

where  C  is  mean  total  cloud  cover  during  daytime  observations 
in  octa. 

2.2.2.  Model  2:  Paltridge  model 

Paltridge  and  Proctor  suggested  a  new  model  that  is  able  to 
determine  the  instant  and  total  daily  radiation  at  any  location  of 
interest.  This  model  takes  the  solar  zenith  angle  (0),  day  length 
(S0)  and  cloud  factor  (CF)  as  inputs.  Their  model  assumes  that  the 
effect  of  atmospheric  water  vapor,  regional  albedo  and  aerosol 
optical  air  mass  on  surface  radiation  is  small  (less  than  5%)  [72]. 

The  hourly  beam  (2b)  and  diffuse  (/d)  radiation  in  (MJ/m2  h)  are 
determined  by: 

Ib  =  3.42286[l—exp(— 0.075(90— 0))]  (131a) 

I  d  =  0.00913  +  0.0125(90— 0)  +  O.723CF  (131b) 

The  mean  monthly  Total  daily  global  radiation  on  horizontal 
surface  (MJ/m2  day)  were  computed  from  correlations  below: 

H  =  Hb+Hd  (131c) 

/•sunset  /•  sunset 

H  =  (1— CF)  /  Ib(6)cosddt+  /  Id(8)dt  (131d) 

J  sunrise  ./sunrise 

In  this  work,  the  integral  was  simply  converted  to  summation 
and  added  up  every  15  min  to  obtain  total  daily  global  radiation. 

The  cloud  factor  (CF)  varies  from  zero  for  clear  sky  to  1  for 
overcast  sky.  This  parameter  could  be  obtained  by  use  of  the 
numbers  of  cloudy  days  in  each  month  and  the  cloud  cover.  Cloud 
cover  is  observed  every  three  hours  in  three  different  ranges:  (0-2) 
octas,  (3-6)  octas,  and  (7-8)  octas.  To  convert  the  cloud  cover  to 


cloud  factor  (CF),  the  following  relationship  is  used  [73]: 

CF  =  (n1)+4.5(n2)  +  7.5(n3) 

8(n3  +n2  +  n3) 


(131e) 


where  n , ,  n2  and  n3  are  the  total  number  of  days  in  each  month, 
with  zero  to  2/8,  3/8  to  6/8,  and  7/8  to  8/8  octas,  respectively. 


2.2.3.  Model  3:  Daneshyar  model 

Following  Paltridge  and  Proctor’s  [72]  work,  Daneshyar  pro¬ 
posed  his  method  by  defining  new  coefficients  for  diffuse  radia¬ 
tion  adjusted  for  the  climate  conditions  of  Tehran,  Iran  [74]: 

The  hourly  beam  and  diffuse  radiation  is  calculated  from 
following  equations: 

Ib  =  3.42286[l—exp(— 0.075(90— 0))]  (132a) 

id  =  0.005 15  +  0.00758(90—  0)  +  0.43677CF  (132b) 

The  mean  monthly  total  radiation  in  (MJ/m2  day)  estimated, 
using  time  steps  of  15  min,  is: 

/•sunset  /*  sunset 

H  =  (1— CF)  /  Ib(8)cos8dt+  /  Id(8)dt  (132c) 

Jsunrise  ,/sunrise 


2.2.4.  Model  4:  Badescu  model 

Badescu  proposed  the  following  correlations  [75]: 

rr-  =  a+bC  (133a) 

no 

4-  =a+bC+cC2  (133b) 

Ho 

=a+bC+cC2  +  dC3  (133c) 

Ho 

2.2.5.  Model  5:  Modified  Daneshyar  model  (Sabziparvar  model) 
Sabziparvar  [76]  in  2007  made  the  following  modifications  to 

the  Daneshyar  method  [74]: 


(1)  In  Daneshyar  method,  solar  constant  of  1353  (W/m2)  has 
been  used  by  the  workers.  Since  the  new  suggested  average 
value  of  solar  constant  is  about  1367  (W/m2)  [77],  the  total 
daily  global  radiation  was  multiplied  by  a  factor. 

(2)  For  each  month,  the  monthly  mean  global  radiation  was 
modified  by  the  Sun-Earth  distance  correction  factor. 

(3)  Height  effect  was  also  applied  separately  on  beam  and  diffuse 
radiation.  For  calculation  of  the  height  correction  factor  in  the 
new  method,  Tehran  (where  Daneshyar  calibrated  his  sug¬ 
gested  model)  is  taken  as  the  reference.  Details  of  this  method 
are  described  in  [76]. 


2.2.6.  Model  6:  Modified  Paltridge  model  (Sabziparvar  model) 

Similar  to  method  5,  Sabziparvar  [76]  in  2007  made  the  same 
modifications  to  the  Paltridge-Proctor  method  [72].  For  locations 
where  sunshine  duration  are  not  observed,  prediction  by  mod¬ 
ified  Paltridge  model  which  only  requires  cloud  data  (are  avail¬ 
able  easily  by  satellites  and  ground-based  measurements)  can  be 
a  good  alternative.  Details  of  this  method  are  described  in  [76]. 

2.3.  Temperature-based  models 

Cloud  observations  and  sunshine  data  are  not  readily  available 
in  all  of  the  locations.  Therefore,  developing  some  precise  solar 
radiation  models  which  use  commonly  available  measured  para¬ 
meters  such  as  air  temperature  is  necessary.  Due  to  the  common 
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availability  of  daily  maximum  and  minimum  air  temperatures, 
several  empirical  methods  have  been  proposed  to  estimate  solar 
radiation  from  these  variables,  especially  for  locations  where  the 
air  daily  temperature  is  the  only  available  meteorological  data. 
The  temperature  based  models  assume  that  the  difference  in 
maximum  and  minimum  temperature  is  directly  related  to  the 
fraction  of  extraterrestrial  radiation  received  at  the  ground  level. 
However,  there  are  factors  other  than  solar  radiation  that  can 
influence  the  temperature  difference  such  as  cloudiness,  humid¬ 
ity,  latitude,  elevation,  topography,  or  proximity  to  a  large  body  of 
water  [78].  In  the  following  section,  the  radiation  models  which 
use  the  maximum  and  minimum  air  temperatures  are  presented 
and  classified  based  on  their  years  of  appearance. 


2.3.1.  Model  1;  Hargreaves  model 

Hargreaves  and  Samani  recommended  a  simple  equation  to 
estimate  solar  radiation  using  only  maximum  and  minimum 
temperatures  [79]: 

tj-  =  a(Tmax—Tmin)05  (134) 

Ho 

Initially,  coefficient  a  was  set  to  0.17  for  arid  and  semi-arid 
regions.  Hargreaves  [80]  later  recommended  using  a  =  0.16  for 
interior  regions  and  a=0.19  for  coastal  regions.  One  of  the 
implications  of  this  model  for  estimating  daily  global  solar 
radiation  is: 


-  Bayat  and  Mirlatifi  model  for  Shiraz,  Iran  [81]: 

7T~  =  0.16(Tmax— Tmjn)0'5  (135) 


2.3.4.  Model  4:  Donatelli  and  Campbell  model 

Donatelli  and  Campbell  algorithm  is  similar  in  structure  to  the 
Bristow  and  Campbell  model  [82],  except  taking  into  account  two 
corrective  functions  of  mean  and  minimum  air  temperature.  The 
model  has  the  following  forms  [84]: 

(a)  ^  =a[l-exp^-h/(ral,g)Ar2/(rmin))]  (139a) 

Tavg  =  (rmax  +  rmin)  (139b) 

f{Tav g)  =  0.017exp(exp(-0.053Ta„g))  (139c) 


f(Tmm)  =  exp(Tmin/c) 


(*»  =  a 

Ho 


1  —exp  I  —b 


ATC\ 
'AT nj 


where  ATm  is  monthly  mean  A T(°C) 
(c)  =  a[l-exp(-fa/(Tal,g)ATc)] 


f  (Tavg)  =  0.017exp(exp(-0.053Tm,gATc)) 


(139d) 

(140) 


(141a) 

(141b) 


2.3.5.  Model  5:  Hunt  et  al.  model 

Hunt  et  al.  [85]  proposed  following  model  by  adding  another 
coefficient  (fa)  to  Hargreaves  and  Samani  model  [79]. 

H  =  o(Tmax—Tmjn)05H0  +  b  (142) 


2.3.2.  Model  2:  Bristow  and  Campbell  model 

Bristow  and  Campbell  developed  a  simple  model  for  daily 
global  solar  radiation  with  a  different  structure  in  which  H  is  an 
exponential  function  of  AT  [82]: 

^  =a[l-exp(-faATc)]  (136) 

where  AT  is  the  temperature  term  difference.  Although 
coefficients  a,  fa  and  c  are  empirical,  they  have  some  physical 
meaning.  Coefficient  a  represents  the  maximum  radiation  that 
can  be  expected  on  a  clear  day.  Coefficients  fa  and  c  control  the 
rate  at  which  a  is  approached  as  the  temperature  difference 
increases. 


2.3.3.  Model  3:  Allen  model 

Following  the  work  of  Hargreaves  and  Samani  [79]  Allen  in 
1997  suggested  employing  a  self-calibrating  model  to  estimate 
mean  monthly  global  solar  radiation  [78]: 

^-=^r(Tmax-Tmin)0-5  (137) 

H  o 

Previously,  Allen  [83]  had  expressed  the  empirical  coefficient 
(I<r)  as  a  function  of  the  ratio  of  atmospheric  pressure  at  the  site 
(Ps,  kPa)  and  at  sea  level  (P0,  101.3  kPa)  as  follows: 


For  the  empirical  coefficient  km,  Allen  suggested  values  of  0.17 
for  interior  regions  and  0.20  for  coastal  regions.  Allen  reported 
that  Eq.  (138)  performs  poorly  for  sites  having  an  elevation  of 
more  than  1500  m. 


2.3.6.  Model  6:  Goodin  et  al.  model 

Goodin  et  al.  refined  the  Bristow  and  Campbell  model  [82]  by 
adding  an  H0-term  meant  to  act  as  a  scaling  factor  allowing  AT  to 
accommodate  a  greater  range  of  H  values  [86]: 


H 

Ho 


=  a 


1-exp 


043) 


The  results  proved  that  this  model  provides  reasonably  accu¬ 
rate  estimates  of  irradiance  at  non-instrumented  sites  and  that 
the  model  can  successfully  be  used  at  sites  away  from  the 
calibration  site  [87], 


2.3.7.  Model  7:  Thorton  and  Running  model 

The  method  proposed  by  Thorton  and  Running  is  based  on  the 
Bristow  and  Campbell  [82]  study  and  is  as  follows  [88]: 

TT -  —  Tr,  ma x?f,  max  (144) 

Ho 

where  it,  max  is  the  maximum  (cloud-free)  daily  total  transmit¬ 
tance  at  a  location  with  a  given  elevation  and  depends  on  the 
near-surface  water-vapor  pressure  on  a  given  day  of  the  year,  and 
t y,  max  stands  for  the  proportion  of  rt,  max  observed  on  a  given  day 
(cloud  correction). 


2.3.8.  Model  8:  Meza  and  Varas  model 

Meza  and  Varas  assumed  that  a  and  c  coefficients  of  Bristow- 
Campbell  [82]  model  are  fixed  and  the  only  fa  coefficient  was 
adjusted  to  minimize  the  square  errors  [89]: 

4.  =0.75[l-exp(-bAr2)] 


(145) 
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2.3.9.  Model  9:  Winslow  model 

The  method  proposed  by  Winslow  et  al.  was  designated  as  a 
globally  applicable  model  and  the  prediction  equation  is  [90]: 


H 

Ho 


=  t  Dl 


1  -a 


^s(Tmin) 


^s(Tmax). 


(146) 


where  es(Tmin)  and  es(Tmax)  are  saturation  vapor  pressures  at  min 
and  max  temperature,  respectively.  Variable  t  accounts  for  atmo¬ 
spheric  transmittance  and  is  estimated  from  site  latitude,  eleva¬ 
tion  and  mean  annual  temperature.  Function  Dl  corrects  the  effect 
of  site  differences  in  day  length,  which  causes  a  variation  between 
the  time  of  maximum  temperature  (and  minimum  humidity)  and 
sunset  [91  ]. 

The  day-length  correction  (Dl)  is  approximated  by: 


Dl  = 


1  Haay 


(ft/4)2 

2  Hiay2 


-1 


(147) 


—  —  a  ln(Tmax— Tmin)+b  (153b) 

Ho 


2.3.14.  Model  14:  Abraha  and  Savage  model 

Abraha  and  Savage  fixed  a=0.75  and  c=  2  in  Donatelli  and 
Campbell  [84]  model  and  proposed  the  following  correlation  [96]: 


£-0-75 


1-exp 


(154) 


2.3.15.  Model  15:  Abraha  and  Savage,  Weiss  et  al.  model 

Weiss  et  al.  and  Abraha  and  Savage  put  a= 0.75  in  Donatelli 
and  Campbell  model  [84]  and  reported  the  following  relationship 
[92,96]: 

^  =0.75[l-exp(-b/(Tal,g)AT2/(Tmin))]  (155) 


2.3.10.  Model  10:  Weiss  et  al.  model 

Weiss  et  al.  simplified  the  Donatelli  and  Campbell  [84]  and 
Goodin  [86]  models  which  only  needs  one  parameter  ( b )  to  be 
calibrated,  while  other  coefficients  (a  and  c)  were  fixed  as  a= 0.75 
and  c=2  [92]: 

(a)  ^=0.75[l-exp(-b/(Tavg)AT2)]  (148a) 

f(Tav g)  =  0.01 7exp  (exp  (-0.053Tm,gAT2))  (148b) 


2.3. 1 6.  Model  1 6:  Almorox  et  al.  model 

A  new  method  is  proposed  by  Almorox  et  al.  for  estimating 
daily  global  solar  irradiance  in  the  form  of  [87]: 

^  =a(Tmax-Tmin),,[l-exp(-c(es(Tmin)/es(Tmax))d)]  (156) 

where  es(Tmin)  and  es(Tmax)  are  saturation  vapor  pressures  at 
minimum  and  maximum  temperature,  respectively.  Saturation 
vapor  pressure  is  related  to  air  temperature  and  the  relationship 
is  expressed  by  [97]: 

es(T)  =  0.61 08exp[17.27(T)/(T+ 273.3)]  (157) 


(b) 


H 

Ho 


0.75 


1-exp 


(149) 


2.4.  Other  meteorological  parameters-based  models 


2.3.11.  Model  11:  Annandale  model 

Annandale  et  al.  modified  Hargreaves  and  Samani  [79]  model  by 
introducing  a  correction  factor  for  parameter  a  to  account  the 
effects  of  reduced  altitude  and  atmospheric  thickness  on  H  as  [93]: 

^  =  fl(l  +2.7  x  10-5z)(Tmax-Tmin)05  (150) 

where  Z  is  the  elevation  above  sea  level  in  m. 

Bayat  and  Mirlatifi  used  this  model  with  a  =  0.1 5  to  estimate 
the  daily  global  solar  radiation  in  Shiraz,  Iran  [81]: 

^  =0.15(1+2.7  x  10-5z)(Tmax-Tmin)05  (151) 


Accurate  prediction  of  the  actual  value  of  solar  radiation  for  a 
given  location  requires  long-term  average  meteorological  data 
which  are  still  scarce  especially  for  underdeveloped  and  devel¬ 
oping  countries.  Therefore  it  is  not  always  possible  to  predict  the 
solar  irradiance  for  a  particular  location  of  interest.  Many 
researchers  have  tried  to  use  various  available  meteorological 
parameters  such  as  precipitation,  relative  humidity,  dew  point 
temperature,  soil  temperature,  evaporation  and  pressure  along¬ 
side  with  the  classical  estimators  such  as  sunshine,  air  tempera¬ 
ture  and  cloudiness  to  predict  the  amount  of  global  solar 
radiation.  In  the  following  section,  the  correlations  which  use 
the  various  meteorological  variables  are  reported  and  classified 
based  on  their  developing  year. 


2.3.12.  Model  12:  Mahmood  and  Hubbard  model 

Mahmood  and  Hubbard  estimated  daily  solar  radiation  based 
on  maximum  and  minimum  daily  air  temperatures  and  proposed 
the  following  model  [94]: 

H  =  a(Tmax— Tmin)0  69H00'91  (152a) 


_  H- 2.4999 
mod  ~  0.8023 


(152b) 


where  Hmod  is  estimated  solar  radiation  corrected  for  systematic 
bias,  in  MJ/m2  day. 


2.3.13.  Model  13:  Chen  et  al.  model 

Chen  et  al.  presented  the  following  models  [95]: 


H 

Ho 


=  a(Tmzx-Tmin)0-5  +  b 


(153a) 


2.4.1.  Model  1:  Swartman  and  Ogunlade  model 

Swartman  and  Ogunlade  stated  that  the  global  radiation  can 
be  expressed  as  a  function  of  the  (S/S0)  ratio  and  mean  relative 
humidity  ( RH )  [98]: 

H  =  a(J^j  RHC  (158a) 

4-  =a  +  b(f)+cRH  (158b) 

Ho  \SoJ 


2.4.2.  Model  2:  Sabbagh  model 

Sabbagh  proposed  following  equation  to  estimate  the  monthly 
average  daily  global  radiation  that  might  be  applicable  to  dry  arid 
and  semi-arid  regions  [99]: 


H  =  0.06407  (Kg')  exp 


\l(s 

RH°  333 

-  1  ^1 

l  \12 

100 

TVnax  ) 

(159) 
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where  I<g  is  geographical  and  seasonal  factor  and  can  be  calcu¬ 
lated  from  the  following  equation: 

Kg  =  100(l ^USo  +  ^m™si)  (160) 

where  is  a  seasonal  factor  in  month  that  is  suggested  by  Reddy 
in  [100]. 

The  model  is  only  reliable  for  locations  with  an  average  mean 
sea  level  of  about  300  m  and  needs  to  be  modified  for  arid  regions 
with  higher  altitudes  [99], 

2.4.3.  Model  3:  Gariepy’s  model 

Gariepy  has  reported  that  the  empirical  coefficients  a  and  b  in 
the  Angstrbm-Prescott  [4]  model  are  dependent  on  mean  air 
temperature  (T,  °C)  and  the  amount  of  mean  precipitation  (P,  cm) 
[101]: 

a  =  0.3791  -0.0041  T-0.01 76  P  (161a) 

b  =  0.4810  +  0.0043  T+0.0097P  (161b) 


2.4.4.  Model  4:  Garg  and  Garg  model 

Garg  and  Garg  proposed  a  double  linear  relation  for  obtaining 
monthly  mean  daily  global  solar  radiation  as  follows  [102]: 

=a+b(-^)+cw  (162a) 

n0  V>o/ 


where  w(cm)  is  the  atmospheric  precipitable  water  vapor  per  unit 
volume  of  air  and  is  computed  according  to  Leckner  [103]: 


w=  0.0049RH 


'exp(26.23— 5416/Tk)’ 

.  T«  . 


(162b) 


2.4.8.  Model  8:  De  Jong  and  Stewart  model 

De  Jong  and  Stewart  introduced  the  effect  of  precipitation  in  a 
multiplicative  form  as  follow  [106]: 

^=a(Tmax-Tmin)b(l+cP+dP2)  (167) 

where  P  is  precipitation  in  mm. 

2.4.9.  Model  9:  Abdalla  model 

Abdalla  modified  the  Gopinathan  [49]  model  for  Bahrain  as  [107]: 
£-=a+b(f)+cT+dRH  (168a) 

no  V>0/ 

=a+b(f]+cT+dRH+ePS  (168b) 

no  VW 

Maghrabi  [108]  estimated  the  ability  of  this  model  for  esti¬ 
mating  monthly  mean  global  solar  radiation  in  kW  h/m2  for 
Tabouk,  Saudi  Arabia: 

=  —0.1 07 +  0.70^  0.0025T+0.004RH  (169) 

no  V>0/ 


2 A A0.  Model  10:  Ododo  et  al.  model 

Ododo  et  al.  proposed  two  new  models  as  follow  [109]: 

b 


tt-  =  a(  Jr- )  WRH' 


Ho 


So 


H 

Ho=e+f 


+gT  max  +  bRH  4-  iTn 


(170a) 

(170b) 


-  El-Metwally  model  for  Egypt  [7] 

=0.219  +  0.526^  +0.004  w  (163) 

no  \So  / 


2.4.1 1.  Model  11:  Hunt  et  al.  model 

Hunt  et  al.  introduced  precipitation  (P,  mm)  and  (Tmax>  °C)  in 
an  additive  form  that  have  had  the  highest  accuracy  at  eight  sites 
in  Canada  [85]: 

H  =  a(rmax— Tmin)  H0  +  bTmix  +  cP+dP  4-  e  (171) 


2.4.5.  Model  5:  Lewis  model 

Lewis  reported  that  global  solar  radiation  on  a  horizontal 
surface  can  be  calculated  by  the  following  equations  [104]: 

H  =  aRHb:  logH  =  a+b  logRH  (164a) 


2.4.12.  Model  12:  Supit  and  Van  Kappel  model 

Supit  and  Van  Kappel  proposed  a  simple  method  to  estimate 
daily  global  radiation  and  tested  it  for  various  locations  in  Europe, 
ranging  from  Finland  to  Italy  [70]. 

H  =  Ho[0\/(Tmax— Tmjn)  +  b\j  (1— C/8)]  +  c  (172) 


H  =  aSbRHc ;  logH  =  a  +  b  logS + c  logRH  ( 1 64b) 

H  =  aexp(bS );  lnH  =  a  +  bS  (164c) 

H  =  a  exp(bRH);  lnH  =  a+bRH  (164d) 

H  =  a  expb(S-RH):  lnH  =  a+b(S—RH )  (164e) 


2.4.6.  Model  6:  Ojosu  and  Komolafe  model 

Ojosu  and  Komolafe  proposed  the  following  equation  [105]: 


H 

TT0=a+ 


b 


+  d 


(165) 


2.4.7.  Model  7:  Gopinathan  model 

Gopinathan  introduced  a  multiple  linear  regression  equation 
of  the  form  [49]: 

=a+bcos(p  +  cZ+d(§-)+eT+fRH  (166) 

Ho  \So  / 


-  Supit  and  Van  Kappel  model  for  London  weather  station,  U.K.  [70] 
H  =  Hot  0.061  ^/(Tmax—Tmin)  +  0.477  (1  C/8)]— 0.557  (173) 


The  regression  coefficients  of  proposed  method  for  various 
locations  in  Europe  are  presented  in  [70]. 

2.4.13.  Model  13:  Togrul  and  Onat  model 

Togrul  and  Onat  investigated  the  effect  of  geographical, 
meteorological  and  astronomical  parameters  on  the  monthly 
mean  global  solar  radiation  and  gave  the  following  correlations 
[110]: 

H  =  a+b  sind  (174a) 

H  =  a+bH0  (174b) 

H  =  a  +  bH0+c(j£)  (174c) 
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H  =  a  +  b  +  c  sinS 


H  =  a  +  b  +  c  sin<5  4-  dT 

H  =  a  +  bH0  +  c(jPj+dST 
H  =  a+bH0+c(jp\  +  dRH+eST 
H  =  a  +  b  sinS+c  )  +  dRH + eT 


(174d) 

(174e) 

(174f) 

(174g) 

(174h) 

(174i) 

(174j) 

(1741<) 

(175a) 

(175b) 

(175c) 

(175d) 


H  =  2.6484  +  2.034^—  )  +  5.2184  sin<5-0.0246T  (175e) 


H  =  — 1.489  +  0. 528H0  +  2. 537  — 0.0067ST  (1750 


H  =  — 0.1764  +  0. 523H0  +  2. 031  "  -0.0146PH-0.0162ST 


(!)■ 

H  =  a  +  bH0  +  c  ^  j  +  d  sinS  +  eRH  +fT 

H  =  a  +  bH0  +  c(J^j+  dRH+eST+fT 

H  =  a  +  bH0  +  c  (JPj  +dsin8  +  eRH+fST+gT 

-  Togrul  and  Onat  model  for  Elazig,  Turkey  [110] 
H  =  4.0891  +  6.459  sinS 


H  =  — 1.32  +  0. 6757H0 
H  =  — 1.3876  +  0. 518H0  +  2.3064 


H  =  2.765  +  2.2984^— J  +4.9597  sinS 


So. 


(175g) 


H  =  3.895  +  5.1 16sin  <5  +  2.44|j+  )  — 0.0143RH— 0.032T 


(175h) 


H  =  5.25— 0.1 69Ho  + 2.501 7^^-^  +6.525  sinS— 0.01 44RH—0.034T 


H  =  — 0.632 +  0.51 6H0  + 2.3255  — 0.0103RH 

+  0.0373ST— 0.062T 


H  =  4.591— 0.1 135H0  +  2.522  +6.1589  sinS 

— 0.01 24i?H  +  0.01 87ST— 0.052T 


(1751) 


(175j) 


(175k) 


2.4.14.  Model  14:  Ertekin  and  Yaldiz  model 

Ertekin  and  Yaldiz  estimated  the  monthly  average  daily  global 
radiation  by  some  multiple  linear  regression  models  using  nine 
different  variables:  extra  terrestrial  radiation,  solar  declination, 
mean  relative  humidity,  ratio  of  sunshine  duration,  mean  tempera¬ 
ture,  mean  soil  temperature,  mean  cloudiness,  mean  precipitation 


(cm)  and  mean  evaporation  (cm)  [2]: 

H  =  a+bS  (176a) 

H  =  a+bHo+cT  (176b) 

H  =  a  +  bH0  +  cf^j+dC  (176c) 

H  =  a+bH0  +  cRH+d(Jpj+eC  (176d) 

H  =  a+bH0  +  cRH+d(J^j+eT+fC  (176e) 

H  =  a+bH0  +  cRH+d(jPj+eT+fC+iE  (176f) 

H  =  a  +  bH0  +  cd  +  dRH+e  +fST+gC+  hE  (176g) 

H  =  a+bH0  +  cd+dRH+e(Jpj+fT+gST+hC+iE  (176h) 

H  =  a  +  bH0  +  cd  +  dRH  +  e(^j+fT+gST+hC+iP+jE  (176i) 

-  Ertekin  and  Yaldiz  model  for  Antalya,  Turkey  [2] 

H=  13.58  +  0.333(5  (177a) 

H  =  — 4.46  +  0.477Ho  +  0.226T  (177b) 

H  =  -2.54  +  0.491Ho  +  4.99  (Jj  — 0.406C  (177c) 


H=  -6.158  +  0.487H0  +  0.0772RH+4.508^|-j-0.633C  (177d) 

H  =  — 7.248  +  0.506H0  +  0.1 13PH+5.987 

— 0.075T— 0.9204C  (177e) 

H  =  — 16.1 6  + 0.493H0  + 0.21 9RH  + 9.282 

— 0.247T — 0.831  C  +  0.2273E  (177f) 

H  =  — 11.90  +  0.353H0  +  0.109<5  +  0.229RH 

+ 11.095  (4^-0.279ST-0.995C+0.241£  (177g) 

H  =  — 12. 57 +  0.361  H0 +  0.0997<5 +  0.234RH+ 1 1.19 

— 0.069T— 0.227ST— 0.99C  +  0.251E  (177h) 

H  =  - 1 3.08  +  0.386H0  +  0.0902  +  0.2254 RH+ 1 1 .59  (Jpj 

—0.034T  +  0.25 1ST— 0.977C—0.0072P+0.2373E  (177i) 

Menges  et  al.  calibrated  Eq.  (176i)  for  Konya,  Turkey  and 
obtained  the  empirical  coefficients  of  this  model  as  follow  [111]: 

H  =  20.296019— 0.096134H0  + 0.3 17593d 

— 0.1 46422RH+ 10.705159  (J^j 

— 0.288332T+0.021331ST+0.359791C 
+  0.207588P— 0.076444E  (178) 
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2.4.15.  Model  15:  Trabea  and  Shaltout  model 

Trabea  and  Shaltout  introduced  the  following  correlation  to  cal¬ 
culate  the  daily  global  solar  radiation  at  five  stations  in  Egypt  [112]: 

±L=a+b(£\+cTmax+dV+eRH+fPS  (179) 

-  Trabea  and  Shaltout  model  for  Egypt  [112] 

^  =  -0.139  +  0.229^-^  +0.009Tmax  + 0.004  1/ 

+  0.002RH  +  0.002PS  (180) 


fraction  (H/H0)  and  various  meteorological  parameters  [19]: 

+  cT  (186a) 


H 

w,-*b 


&■ 


H 

Hi 


+  cRH 


H 

Ho 

H 

Ho 


=  a+bT+cRH 
=  u+b(Tmax— Tmin)  +  cC 


(186b) 

(186c) 

(186d) 


2.4. 3  6.  Model  3  6;  El-Metwally  model 

El-Metwally  proposed  following  simple  new  methods  to  esti¬ 
mate  global  solar  radiation  based  on  meteorological  data  over  six 
stations  in  Egypt  [113]: 


H  —  aH0  +  frTmax  +  cTmin  +  c/C+ e 

(181a) 

H  =  qHo  +  bTmax-\-  cTm\n-\-  dC 

(181b) 

H  =  exp(  aH0  +  bT  max  +  cTmin  +  dC  +  e) 

(181c) 

2.4. 3  7.  Model  3  7;  Chen  et  al.  model 

Chen  et  al.  suggested  a  logarithmic  relationship  between  the 
daily  global  solar  radiation,  daily  extraterrestrial  solar  radiation 
and  the  temperature  difference  between  the  maximum  and 
minimum  daily  air  temperature  as  follow  [95]: 

—  =  a  ln(Tmax—Tmjn)+b^-J  +d  (182) 

2.4.18.  Model  18:  Modified  Sabbagh  method  (Sabziparvar  model) 
Sabziparvar  in  2007  revised  Sabbagh  [99]  model  to  predict  the 

monthly  average  daily  solar  radiation  on  horizontal  surfaces  in 
various  cities  in  central  arid  deserts  of  Iran.  The  modification  was 
made  by  the  inclusion  of  following  factors  to  Sabbagh  model: 

(1)  Height  correction  factor,  (2)  Sun-Earth  distance  correction 
factor  and  (3)  inclusion  of  monthly  total  number  of  dusty  days  [76], 

2.4.19.  Model  19:  Bulut  and  Biiyilkalaca  model 

Bulut  and  Buyiikalaca  proposed  a  simple  model  for  estimating 
the  daily  global  radiation  on  a  horizontal  surface.  The  model  is 
based  on  a  trigonometric  function,  which  has  only  one  indepen¬ 
dent  parameter,  namely  the  day  of  the  year  [114]: 

1=  b+(a-b)|sin[=^=(n  +  5)]  |  (183) 

They  examined  their  model  for  68  provinces  of  Turkey  [114], 
Some  investigator  employed  this  model  to  estimate  global 
solar  radiation  in  different  locations  as  follows: 


—  o+b(Tmax— Tmjn)05  +  cC  (186e) 

m-°+b(k)+*  1,860 

-  El-Sebaii  et  al.  model  For  Jeddah,  Saudi  Arabia  [19] 

^  =  -1.92  +  2.60K+) +  0.006 T  (187a) 

no  \-->0  / 

^  =  -1.62  +  2.24^+0.332 RH  (187b) 

no  V-W 

=0.139-0.0031+  0.896RH  (187c) 

no 

£-  =0.214  +  0.035(Tmax-rmin)-0.028C  (187d) 

no 

=  — 0.08 +  0.2  l(Tmax—Tmin)05— 0.01 2C  (187e) 

no 

^-  =  -2.76  +  3.72^)+0.001  C  (187f) 

Ho  V+o/ 

2.4.21.  Model  21:  Maghrabi  model 


Maghrabi  established  a  simple  model  to  calculate  the  monthly 
mean  global  solar  radiation  on  a  horizontal  surface  using  five 
meteorological  parameters  [108]: 

H  =  a+b(Jp\+cT+dP0+ePWV+fRH  (188) 

The  precipitable  water  vapor  ( PWV)  is  calculated  using  the 
Reitan  [115]  equation  as  follows: 

In  PWV  (mm)  =  0.1102  +  0.0613  td  (189) 

-  Maghrabi  model  for  Tabouk,  Saudi  Arabia  [108] 

H(kW  h/m2)  =  163.01-1. 04 (jL)  +0.12T-0.21P0 

--l.06PWV-0.03RH  (190) 


-  Bulut  and  Btiytikalaca  model  for  Ankara,  Turkey  [114] 

3  =  3.86  +  (22.71-3.86)|sin[^(n  +  5)]|15  (184) 

-  El-Sebaii  et  al.  model  For  Jeddah,  Saudi  Arabia  [19] 

3  =  14.92  +  9.61  |sin  [3JJ5  (n+5)]  (' *  (185) 

2.4.20.  Model  20:  El-Sebaii  et  al.  model 

El-Sebaii  et  al.  developed  and  presented  empirical  correlations 
between  the  monthly  average  of  daily  global  solar  radiation 


3.  A  case  study 

3.3.  Experimental  data 

It  is  already  mentioned  that  the  main  objective  of  this  study  is 
to  comprehensively  collect  and  review  the  global  solar  radiation 
models  available  in  the  literature  and  categorize  them  based  on 
the  employed  meteorological  parameters.  In  order  to  evaluate  the 
applicability  and  accuracy  of  the  collected  models  for  computing 
the  monthly  average  daily  global  solar  radiation  on  a  horizontal 
surface,  the  geographical  and  meteorological  data  of  Yazd,  Iran 
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was  used.  Yazd  city  is  located  in  latitude  31‘54'N,  longitude 
54°17'E  and  its  altitude  of  sea  level  is  1237.2  m.  Yazd  province 
is  situated  in  the  central  part  of  Iranian  plateau  and  because  of 
being  in  the  sun-belt,  it  is  an  ideal  location  to  benefit  the 
advantages  of  solar  energy  utilization  and  adoption  of  its  related 
technologies  [116].  The  total  radiation  received  by  a  horizontal 
surface  in  Yazd  City  during  a  year  is  reported  to  be  around 
7787  MJ/m2  with  about  3270  sunshine  hours  in  a  year  while  the 
number  of  overcast  days  is  fewer  than  1110  h  per  year  [117].  In 
the  present  work,  the  value  of  measured  data  of  daily  global  solar 
radiation  on  a  horizontal  surface  in  Yazd  city  in  the  period  of 
1982-2008,  as  well  as  number  of  sunshine  hours,  maximum, 
minimum  and  daily  ambient  temperatures,  amount  of  cloud 
cover,  relative  humidity  and  the  other  related  meteorological 
parameters  during  the  period  of  1988-2008  were  obtained  from 
the  Islamic  Republic  of  Iran  Meteorological  Office  (1RIMO)  data 
centre. 

The  measured  global  solar  radiation  data  were  checked  and 
controlled  for  errors  and  inconsistencies.  The  purpose  of  data 
quality  control  was  to  eliminate  faulty  data  and  inaccurate 
measurements.  After  the  quality  control,  the  measured  data 
were  averaged  to  obtain  the  monthly  mean  daily  values  by  taking 
the  data  for  the  average  day  of  each  month  as  recommended 
by  Duffie  and  Beckman  [118],  The  measured  data  were  then 
divided  into  two  sets.  The  first  sub-data  set  (1988-2003)  were 
employed  to  develop  empirical  correlations  between  the  monthly 
average  daily  global  solar  radiation  fraction  (H/H0)  and  monthly 
average  of  desired  meteorological  parameters  (for  the  global  solar 
radiation  this  subset  is  contained  the  measured  data  of  1982- 
2003).  The  second  sub-data  set  (2004-2008)  were  then  used  to 
validate  and  evaluate  the  derived  models  and  correlations. 
The  measured  long-term  monthly  average  daily  global  radiation 
distribution  for  Yazd  city  in  the  periods  of  1982-2008  is  shown 
in  Fig.  1. 


3.2.  Statistical  evaluation 

The  accuracy  and  performance  of  the  derived  correlations  in 
predicting  of  global  solar  radiation  was  evaluated  on  the  basis  of 
the  following  statistical  error  tests  which  are  coefficient  of 
determination  (R2),  root  mean  square  error  ( RMSE ),  mean  bias 
error  ( MBE ),  mean  absolute  bias  error  ( MABE ),  mean  percentage 
error  (MPE),  correlation  coefficient  (r)  and  t-Test  statistic.  These 


error  indices  are  defined  as: 

p2  1 

£"=1  (Hi.m-H^)2 


RMSE  = 


n  £  Mu 


1/2 


MBE 


=  „E 


MABE=lj2  \(HUc-Hlm) 


(191) 

(192) 

(193) 

(194) 

(195) 


E"=i  (Htc-Hc)(H,.m-Hm) 

EL  i  (Hu-H 7)2  EL  i  (Hi.m-7^)2]1'2 


(196) 


In  the  above  relations,  the  subscript  i  refers  to  the  ith  value  of 
the  solar  irradiation  and  n  is  the  number  of  the  solar  irradiation 
data.  The  subscripts  “c"  and  “m”  refer  to  the  calculated  and 
measured  global  solar  irradiation  values,  respectively.  Hc  is  the 
mean  calculated  global  radiation  and  Hm  is  the  mean  measured 
global  radiation. 

Although  these  statistical  indicators  generally  provide  reason¬ 
able  criteria  to  compare  models,  but  do  not  objectively  indicate 
whether  the  estimates  from  a  model  are  statistically  significant. 
The  t-statistic  (Eq.  (197))  allows  models  to  be  compared  and  at 
the  same  time  indicates  whether  a  model’s  estimate  is  statisti¬ 
cally  significant  at  a  particular  confidence  level  or  not. 


t  = 


(n—V)(MBE)2 
(RMSE)2 -(MBE)2 


1/2 


(197) 


The  smaller  the  t-value,  the  better  the  model  performance  is. 
To  determine  whether  the  estimates  from  a  model  are  statistically 
significant,  one  simply  has  to  determine  from  the  standard 
statistical  tables,  the  critical  t-value,  i.e.,  (ta/2)  at  a  level  of 
significance  and  (n-l)°  of  freedom.  The  model  is  judged  to  be 
statistically  significant  if  the  calculated  t-value  is  less  than  the 
critical  value. 


Fig.  1.  The  measured  values  of  monthly  average  daily  global  solar  radiation  for 
Yazd  city  in  the  periods  of  1982-2008. 


3.3.  Result  and  discussion 
3.3.1.  Calibration 

Linear,  multiple  linear  and  nonlinear  regressions  were  carried 
out  between  the  monthly  mean  measured  global  solar  radiation 
and  the  meteorological  parameters  using  the  first  sub-data  set  to 
obtain  the  values  of  empirical  coefficients  of  the  selected  models 
from  each  category.  The  coefficient  of  determination  (R2)  index  is 
used  to  determine  that  how  well  the  regression  line  approximates 
the  real  data  points.  A  model  is  more  efficient  when  R2  is  closer  to 
1.  In  addition,  the  t-value  index  is  used  to  determine  whether  the 
estimates  from  a  model  are  statistically  significant  or  not.  For  the 
model’s  estimates  to  be  judged  statistically  significant  at  the 
(1-a)  confidence  level,  the  calculated  t-value  must  be  less  than 
the  critical  value.  Regression  coefficients  of  the  selected  models  in 
each  category  for  Yazd  city  along  with  the  values  of  coefficients  of 
determination  (R2)  and  t-values  are  presented  in  Table  1.  It  should 
be  mentioned  that  several  radiation  models  from  each  category  is 
used  for  estimating  the  monthly  average  daily  global  solar 
radiation  fraction  (H/H0),  however,  the  results  of  one  selected 
model  from  each  category  is  presented  in  Table  1  for  brevity. 
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Table  1 

Regression  coefficients  of  selected  radiation  models  for  Yazd  city  along  with  the  values  of  coefficients  of  determination  (R2)  and  t-value. 


Models 

(R2) 

t-value 

a 

b 

c 

d 

Sunshine-based  models 

El-Metwally  model  (Eq.  (116)) 

Cloud-based  models 

0.979 

0.5234 

0.7272 

Badescu  model  (Eq.  (133a)) 

Temperature-based  models 

0.9948 

0.2663 

0.7658 

-0.0594 

Hargreaves  model  (Eq.  (134)) 

0.9543 

0.2625 

0.1746 

Other  meteorological  parameters-based  models 

Chen  et  al.  model  (Eq.  (182)) 

0.9971 

0.3782 

-0.01088 

5.127 

0.101 

-4.299 

Table  2 

Statistical  results  for  the  validation  of  the  selected  models  for  Yazd  city  (using  data 
in  the  period  of  2004-2008). 


Models  RMSE 

(MJ/m2 
day) 

MBE 

(MJ/m2 

day) 

MABE 

(MJ/m2 

day) 

MPE 

(%) 

r 

Sunshine-based  models 

El-Metwally  model  0.5385 

0.04371 

0.4124 

0.8485 

0.9969 

Cloud-based  models 

Badescu  model  (a)  1.152 

-0.3385 

0.9592 

-1.877 

0.9846 

Temperature-based  models 

Hargreaves  model  0.7103 

-0.02285 

0.6293 

0.3039 

0.9934 

Other  meteorological  parameters-based  models 

Chen  et  al.  model  0.8542  0.2104  0.746 

0.9859 

0.9913 

Month 


Fig.  2.  Measured  (in  the  periods  of  2004-2008)  and  calculated  values  of  monthly 
average  daily  global  solar  radiation  for  Yazd  city. 


The  results  of  Table  1  show  that  all  the  selected  models  give 
the  coefficients  of  determination  ( R 2)  higher  than  95%,  indicating 
a  very  good  fitting  between  the  monthly  average  daily  global 
radiation  and  the  other  meteorological  parameters.  The  best  fit  is 
obtained  by  employing  Chen  et  al.  model  [95].  Comparison 
between  the  different  models  according  to  the  t-value  shows  that 
the  calculated  t-values  are  less  than  the  critical  t-value  (2.201  at 
95%  confidence  level),  indicating  that  all  models  have  statistical 
significance. 


3.3.2.  Validation 

The  performance  of  developed  correlations  in  predicting  the 
global  solar  radiation  on  horizontal  surface  is  also  assessed  by 
comparing  each  models  outputs  and  the  measured  values  from 
the  second  subset  data  based  on  the  r,  RMSE,  MBE,  MABE  and  MPE 
error  indices.  A  summary  of  all  the  statistical  parameters  are 
presented  in  Table  2.  For  higher  modeling  accuracy  RMSE,  MBE, 


MABE  and  MPE  indices  should  be  closer  to  zero,  but  correlation 
coefficient  (r)  should  approach  to  1  as  closely  as  possible. 

The  performance  of  the  most  accurate  models  in  each  category 
for  estimating  the  monthly  average  daily  global  solar  radiation  in 
Yazd  is  presented  and  compared  in  Table  2.  From  the  statistical 
analysis,  it  can  be  seen  that  the  estimated  values  of  monthly 
mean  daily  global  solar  irradiation  are  in  good  agreement  with 
the  measured  values  for  all  models.  It  is  found  that  the  values  of 
RMSE  are  in  the  range  from  0.5385  to  1.152  (in  MJ/m2  day).  The 
MBE  achieved  in  this  study  for  all  models  are  in  the  acceptable 
range.  Negative  values  of  MBE  in  temperature  and  cloud  based 
models  indicate  an  underestimation  of  measured  global  solar 
radiation  by  these  models.  As  an  over-estimation  of  an  individual 
observation  may  cancel  under-estimation  in  a  separate  observa¬ 
tion,  using  MABE  index  is  more  appropriate  than  MBE.  The  biggest 
value  of  MABE  is  belonged  to  Badescu  model  (a)  with  0.9592.  The 
mean  percentage  errors  (MPE)  of  all  models  are  in  the  range  of 
acceptable  values  between  0.3039%  and  -1.877%.  Also  according 
to  the  statistical  test  of  the  correlation  coefficient  (r),  all  models 
give  very  good  results  (above  0.98).  El-Metwally  sunshine  based 
model  [7]  gives  the  highest  value  of  r  with  value  of  0.9969, 
indicating  a  strong  positive  linear  relationship  between  the 
measured  and  calculated  values  of  the  global  solar  radiation. 

The  measured  values  of  the  monthly  average  daily  global 
radiation  and  corresponding  calculated  values  by  employing 
El-Metwally  model  [7],  Badescu  model  [75],  Hargreaves  model 
[79]  as  well  as  Chen  et  al.  model  [95]  for  Yazd  city  are  illustrated 
in  Fig.  2.  As  may  be  seen,  agreement  between  the  values  obtained 
from  the  selected  models  and  the  measured  data  are  very  good. 

Based  on  the  statistical  indicators  presented  in  Table  1  as  well 
as  Fig.  2,  all  models  show  good  estimation  of  the  monthly  average 
daily  global  solar  radiation  on  a  horizontal  surface  for  Yazd  city. 
But  the  El-Metwally  sunshine  based  model  [7]  produces  the  best 
result  among  the  correlations  developed  in  this  study.  Therefore, 
this  model  is  recommended  for  the  prediction  of  global  solar 
radiation  on  a  horizontal  surface  in  Yazd  and  elsewhere  with 
similar  climatic  conditions,  where  the  radiation  data  are  missing 
or  unavailable. 

The  agreement  between  measured  and  calculated  value  of  H 
was  also  confirmed  for  each  month  by  calculating  the  relative 
percentage  error  (e).  The  relative  percentage  error  for  each  month 
is  defined  as: 

e=  (Hi'[ x  100  (198) 

V  Hi.m  J 

Fig.  3  shows  the  comparisons  of  the  relative  percentage  error 
(e)  of  measured  and  calculated  H  for  all  developed  models.  It  is 
seen  from  Fig.  3  that  the  relative  percentage  error  for  each  month 
rarely  exceeds  +10%. 

It  is  found  from  Fig.  3  that  almost  all  methods  provide  low 
performance  at  both  winter  and  spring,  this  is  due  to  the  increase 
of  cloudiness  and  aerosol  content,  respectively. 
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—♦—El-Metwally  model  [7] 
— Badescu  model  [75] 
—A— Hargreaves  model  [79] 
— Chen  et  al.  model  [95] 


Fig.  3.  Comparisons  of  monthly  relative  percentage  error  (e)  of  monthly  average 
daily  measured  and  calculated  global  solar  radiation  for  Yazd  (2004-2008). 


4.  Concluding  remarks 

Solar  energy  technologies  offer  opportunities  for  use  of  a  clean, 
renewable,  and  domestic  energy  resource,  and  are  essential 
components  of  a  sustainable  energy  future.  Availability  of  a 
complete  and  accurate  solar  radiation  data  base  for  each  specific 
region  is  essential  in  the  prediction,  study,  and  design  of  solar 
energy  systems.  Information  on  global  solar  radiation  at  any  sites 
(preferably  obtained  over  a  long  period  of  time)  is  useful  not  only 
to  the  locality  where  the  radiation  data  is  collected  but  also  for 
the  wider  world  community. 

A  global  study  of  the  world  distribution  of  global  solar 
radiation  requires  knowledge  of  the  radiation  data  in  various 
countries.  For  the  purpose  of  worldwide  marketing,  the  designers 
and  manufactures  of  solar  equipment  will  need  to  know  the 
average  global  solar  radiation  available  in  different  regions.  In  this 
regard,  a  common  practice  is  to  estimate  average  daily  global 
solar  radiation  using  appropriate  correlations  which  are  empiri¬ 
cally  established  by  using  the  measured  data  at  some  selected 
locations.  These  correlations  estimate  the  values  of  global  solar 
radiation  from  more  readily  available  meteorological,  climatolo¬ 
gical,  and  geographical  parameters.  The  number  of  such  correla¬ 
tions  or  empirical  equations  that  have  been  developed  and  tested 
is  relatively  high,  which  makes  it  difficult  to  choose  the  most 
appropriate  correlation  for  a  particular  purpose  and  site. 

In  this  study,  a  comprehensive  collection  and  review  of  available 
solar  models  is  conducted  to  assist  in  the  selection  of  most 
appropriate  and  accurate  model  based  on  the  available  measured 
meteorological  data.  Moreover,  the  collected  models  are  classified 
into  four  categories  namely,  sunshine-based,  cloud-based,  tempera¬ 
ture-based,  and  other  meteorological  parameter-based  models.  The 
regression  constants  of  some  collected  solar  models  have  been 
generally  presented  to  calculate  the  global  solar  radiation  with  high 
accuracy  in  a  given  location.  These  regression  coefficients  have  been 
computed  using  measured  monthly  mean  daily  global  radiation  and 
other  meteorological  variables  for  locations  of  interest.  The  regres¬ 
sion  coefficients  in  the  presented  models  differ  from  location  to 
location.  This  may  be  due  to  the  local  and  seasonal  changes  in  the 
type  and  thickness  of  cloud  cover,  the  effects  of  snow  covered 
surfaces,  the  concentrations  of  pollutants,  and  latitude  [13]. 

The  correlations  reported  in  this  study  have  had  the  best 
estimation  in  the  location  of  interest  and  the  present  classification 
of  solar  models  make  it  easy  to  choose  a  suitable  model  from 
various  existing  models  based  on  the  meteorological  parameters 
measured  in  the  location  of  interest.  These  correlations  with  these 
selected  variables  can  be  used  by  solar  energy  system  designers  in 
the  design  and  performance  prediction  of  solar  applications  with 
fairly  high  level  of  accuracy  in  calibration  sites  and  possibly  in 
regions  with  similar  climatic  conditions. 


Furthermore,  to  evaluate  the  accuracy  and  applicability  of 
collected  various  models  for  estimating  the  monthly  average 
daily  global  radiation  on  a  horizontal  surface,  the  long  term 
measured  meteorological  and  solar  data  of  Yazd  city,  Iran  is 
employed.  With  employing  linear  and  nonlinear  least  square 
method,  regression  coefficients  of  some  selected  models  from 
the  four  categories  were  obtained  by  using  first  sub-data  set.  The 
developed  models  were  then  compared  with  each  other  and  with 
the  experimental  data  of  second  subset  on  the  basis  of  the 
statistical  error  indicators  such  as  RMSE,  MBE,  MABE,  MPE  and 
correlation  coefficient  (r).  The  most  accurate  models  in  each 
category  are  then  identified  and  presented.  These  models  have 
reasonable  values  of  estimation  errors.  Based  on  the  statistical 
results,  El-Metwally  sunshine-based  model  reproduce  the  mea¬ 
sured  monthly  average  daily  global  solar  irradiation  for  Yazd  city 
with  the  highest  accuracy. 

A  common  feature  of  almost  all  models  collected  in  this  paper 
is  that  they  account  for  latitude,  solar  declination,  elevation,  day 
length  and  atmospheric  transmissivity  by  including  the  extra¬ 
terrestrial  radiation  (H0)  term  in  the  model.  Among  the  four 
categories  of  solar  radiation  models,  the  sunshine-based  methods 
are  generally  more  accurate  but  they  are  often  limited  by  the  lack 
of  availability  of  sunshine  records.  Estimation  of  global  solar 
radiation  from  the  air  temperature  offers  an  important  alternative 
in  the  absence  of  measured  H  or  sunshine  duration  because  of  the 
wide  availability  of  air  temperature  data.  The  main  advantage  of 
models  in  this  category  is  the  readily  available  data.  Temperature- 
based  models  are  a  convenient  tool  for  estimating  solar  radiation 
if  the  parameters  can  be  calibrated  for  each  specific  location  [87]. 
The  temperature  range  (AT)  is  the  main  factor  affecting  accuracy 
of  the  temperature-based  models.  Larger  (AT)  generally  results  in 
a  better  predictive  accuracy,  meaning  that  the  temperature-based 
models  are  more  applicable  in  areas  with  larger  temperature 
range.  This  implies  that  model  calibration  is  particularly  sensitive 
in  humid  regions  where  (AT)  is  generally  small  [119]. 

When  sunshine  duration  data  are  not  available,  prediction  by 
the  method  that  only  requires  cloud  data  (are  available  easily  by 
satellites  and  ground-based  measurements)  can  be  a  good  alter¬ 
native.  One  drawback  of  these  models  is  that  they  are  sensible  to 
human  biasing  [70], 

The  models  presented  in  the  last  category,  employ  more  than 
one  meteorological  parameter  for  predicting  the  global  solar 
radiation.  These  models  reported  a  good  estimation  of  solar 
radiation,  but  they  are  limited  because  of  using  the  various 
meteorological  parameters  as  models  inputs  that  are  not  readily 
available  in  most  of  the  location  of  interest.  Moreover,  to  prevent 
the  error  of  measurement  tools  also  increasing  speed,  the  solar 
radiation  should  be  estimated  using  the  minimum  measured 
input  parameters. 
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